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EXECUTIVE  SUMMARY 


Electron  donor  (ED)  addition  can  be  very  effective  in  stimulating  enhanced  reductive 
bioremediation  (ERB)  of  a  wide  variety  of  groundwater  contaminants.  However,  ERB  can  result 
in  Secondary  Water  Quality  Impacts  (SWQIs)  including  decreased  levels  of  dissolved  oxygen 
(O2),  nitrate  (NO3'),  and  sulfate  (SO42"),  and  elevated  levels  of  dissolved  manganese  (Mn2+), 
dissolved  iron  (Fe2+),  methane  (CH4),  sulfide  (S2~),  organic  carbon,  taste  and  odor  concerns,  and 
naturally  occurring  hazardous  compounds  (e.g.,  arsenic  [As]).  Fortunately,  this  ‘plume’  of 
impacted  groundwater  is  usually  confined  within  the  original  contaminant  plume  and  is  unlikely 
to  adversely  impact  potable  water  supplies. 

This  report  presents  the  results  of  research  performed  as  part  of  project  ER-2131  “Numerical 
Modeling  of  Post-Remediation  Impacts  of  Anaerobic  Bioremediation  on  Groundwater  Quality” 
supported  by  the  Strategic  Environmental  Research  and  Development  Program  (SERDP)  and 
performed  by  a  team  of  researchers  from  North  Carolina  State  University  (NCSU)  and  the 
United  States  Geological  Survey  (USGS). 

OBJECTIVES 

The  overall  objective  of  this  research  is  to  develop  an  improved  understanding  of  the  near-  and 
long-term  impacts  to  groundwater  quality  after  implementation  of  in  situ  anaerobic 
bioremediation  processes.  Specific  objectives  achieved  during  this  project  include: 

A.  Formulate  a  general  modeling  approach  appropriate  for  simulating  the  production  and 
subsequent  natural  attenuation  of  SWQIs  for  a  wide  range  of  sites.  The  model  approach  and 
validation  was  based  on  data  from  the  Bemidji  Crude  Oil  Spill  and  Cape  Cod  Wastewater 
Plume. 

B.  Assemble  a  database  of  SWQIs  at  ERB  sites  to  determine  the  range  of  SWQI  concentrations 
and  evaluate  trends  over  time  and  distance  downgradient. 

C.  Use  the  validated  modeling  approach  and  results  from  the  database  analysis  to  develop  a 
general  conceptual  model  of  SWQI  formation,  mobilization,  and  attenuation.  This 
conceptual  model  can  be  used  to  identify  sites  where  SWQI  could  be  an  issue  and  develop 
monitoring  and  management  approaches  to  minimize  adverse  impacts. 

TECHNICAL  APPROACH 

A  reactive  transport  model  was  developed  to  simulate  the  production  and  natural  attenuation  of 
important  electron  acceptors  and  SWQI  parameters  including  organic  carbon,  O2,  NO3',  Mn2+, 
Fe2+,  SO42',  S2",  CH4,  and  CO2.  The  behavior  of  these  and  other  reactive  species  in  groundwater 
have  been  intensively  studied  for  many  years  by  the  USGS  Toxics  Substances  Hydrology 
Program  (Toxics  Program).  Model  simulation  results  were  compared  to  an  extensive  database 
on  monitoring  results  from  two  Toxics  Program  sites:  (a)  Bemidji  Crude  Oil  Spill,  MN;  and  (b) 
Cape  Cod  MMR  Wastewater  Plume,  MA. 
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SWQI  data  were  collected  from  47  active  and  former  ERB  sites  across  the  United  States, 
including  different  contaminants  of  concern  (COCs),  injection  strategies,  organic  substrates,  and 
geographic  locations.  Monitoring  data  was  compiled  on  SWQI  parameters  in  917  different 
monitoring  wells  from  regulatory  monitoring  reports,  DOD  (i.e.,  AFCEE,  ESTCP,  SERDP) 
project  reports,  journal  articles,  and  personal  correspondence  with  site  project  managers. 
Summary  statistics  (mean,  median,  standard  deviation,  etc.)  and  cumulative  frequency 
distribution  plots  were  generated  for  upgradient,  treatment  zone,  and  various  downgradient 
locations  to  define  the  range  of  SWQIs  impacts  and  evaluate  trends  over  time  and  distance 
downgradient. 

The  reactive  transport  model  and  SWQI  database  were  then  used  to  produce  indicator 
simulations  that  illustrate  the  major  processes  control  the  formation  and  natural  attenuation  of 
SWQIs  at  ERB  sites.  The  model  was  used  to  simulate  the  fate  of  the  added  organic  carbon  and 
how  the  reduced  products  affected  the  aquifer  over  a  period  of  40  years  following  subustrate 
addition. 

Results  from  the  reactive  transport  model  development,  SWQI  database,  and  indicator 
simulations  were  integrated  to  develop  a  general  conceptual  model  of  the  major  processes 
controlling  SWQI  production  and  attenuation.  This  conceptual  model  can  be  used  as  a  guide  in 
understanding  the  magnitude,  areal  extent,  and  duration  of  SWQIs  in  ERB  treatment  zones  and 
the  natural  attenuation  of  SWQI  parameters  as  the  dissolved  solutes  migrate  downgradient  with 
ambient  groundwater  flow. 

RESULTS 

Reactive  Transport  Model  Development 

Reactive  transport  models  were  created  for  USGS  long-term  study  sites  where  reducing 
conditions  were  created  in  shallow  aquifers  by  crude  oil  (Bemidji)  or  wastewater  disposal  (Cape 
Cod).  The  goal  of  the  modeling  was  to  understand  the  creation,  natural  attenuation,  and 
longevity  of  the  SWQIs  created  by  the  reducing  conditions.  The  model  simulations  described 
flow  and  reactions  along  two-dimensional  vertical  cross  sections  of  the  study  aquifers.  Steady- 
state  models  of  the  flow  systems  were  created  using  MODFLOW  (Harbaugh,  2005)  and  reactive 
transport  simulations  were  performed  with  PHT3D  (Prommer  et  al.,  2003). 

For  the  Bemidji  Crude  Oil  Spill,  the  goal  of  the  modeling  was  to  understand  the  generation  and 
fate  of  the  most  important  SWQIs  at  the  site,  especially  CEE,  aqueous  Fe(II),  and  dissolved 
organic  carbon  (DOC).  To  constrain  the  model,  complete  datasets  were  assembled  for  four  time 
points  representing  the  time  of  the  oil  spill  in  1979  (i.e.,  immediately  after  the  oil  spill),  1987±2 
years,  1993±2  years,  and  2008±2  years.  The  model  formulation  included  multiple  hydrocarbon 
compounds  in  the  oil,  dissolved  hydrocarbons,  sediment  and  aqueous  Fe  and  Mn,  mineral 
controls  on  pH  and  alkalinity,  and  outgassing  of  CH4  and  carbon  dioxide  (Ng  et  al.,  2015).  The 
simulation  covered  the  time  period  from  1979  to  2008  along  a  260-m  horizontal  by  7-m  vertical 
cross-section.  Other  aspects  of  the  model  including  initial  concentrations,  grid  discretization, 
permeabilites,  dispersivities,  dissolution,  equilibrium  and  kinetic  reactions,  outgassing,  and 
sorption  are  described  in  Ng  et  al.  (2015). 
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For  the  Cape  Cod  Wastewater  Plume,  the  modeling  goal  was  to  understand  what  controlled  the 
recovery  rate  of  aerobic  conditions  in  the  aquifer  after  the  cessation  of  wastewater  disposal.  The 
time  frame  of  recovery  depended  on  re-oxidation  rates  of  two  pools  of  reduced  species  that  had 
accumulated  on  the  aquifer  sediments  during  the  60  years  of  wastewater  disposal  pond 
operations.  These  two  species  are  sediment-bound  biogenic  organic  carbon  located  below  the 
disposal  ponds  and  sorbed  Fe[II]  located  in  the  aquifer  downgradient  of  the  ponds.  The  initial 
geochemical  conditions  in  the  model  were  based  on  observations  at  the  time  of  disposal  cessation 
in  1996.  The  initial  spatial  distribution  of  sediment-bound  organic  carbon  and  reaction  rates 
were  calibrated  to  the  observed  aquifer  recovery  rates  in  the  first  10  years  after  cessation. 
Competitive  sorption  parameters  for  major  cations,  Fe[II],  Mn[II],  and  H+  were  based  on 
laboratory  experiments. 

SWQI  Database 

Analysis  of  the  SWQI  database  results  provided  significant  insights  into  the  production  and 
attenuation  of  SWQI  parameters.  Maximum  TOC  concentrations  decline  rapidly  with  distance 
downgradient  of  the  injection  zone.  As  a  result,  CH4  production  and  SO42'  reduction  are  largely 
restricted  to  the  injection  zone  and  immediately  downgradient. 

Dissolved  Mn  and  Fe  concentrations  increased  in  the  injection  zone  at  essentially  all  sites,  then 
declined  with  distance  downgradient.  At  >50  m  downgradient,  the  median  increase  in  dissolved 
Mn  and  dissolved  Fe  was  0.08  mg/L  and  0.38  mg/L,  respectively.  Arsenic  was  not  monitored  at 
most  sites.  When  monitoring  data  was  available,  arsenic  concentrations  often  increased 
following  ED  addition.  Increasing  dissolved  arsenic  concentrations  were  correlated  with 
increasing  Fe  concentrations  consistent  with  arsenic  release  during  reduction  of  Fe  oxides  and/or 
arsenic  removal  during  sorption/co-precipitation  with  Fe.  In  the  injection  area,  sulfate 
concentrations  were  often  reduced  by  one  to  three  orders  of  magnitude,  but  were  not  completely 
depleted  at  most  sites.  Further  downgradient,  the  change  in  sulfate  concentration  was  smaller, 
indicating  substantial  attenuation  of  the  low  sulfate  plume.  Sulfide  concentrations  typically 
increased  by  less  than  1  mg/L  in  injection  area  wells.  The  maximum  increase  in  CH4 
concentration  declines  with  distance  downgradient,  indicating  some  attenuation.  However,  there 
were  a  significant  number  of  wells  at  greater  than  50  m  downgradient  with  high  CH4  levels  (>10 
mg/L). 

While  SWQIs  were  detected  in  downgradient  wells,  the  large  majority  of  these  wells  were  within 
the  primary  contaminant  plume.  In  wells  where  chlorinated  ethenes  were  previously  below 
maximum  contaminant  levels  (MCLs),  SWQIs  were  generally  within  background  levels  or 
within  a  few  meters  of  the  injection  zone.  As  a  result,  SWQI  production  is  unlikely  to  adversely 
impact  potable  water  supplies. 

Indicator  Simulations 

Results  from  the  indicator  simulations  provided  insights  into  the  major  processes  controlling  the 
formation  and  natural  attenuation  of  SWQIs  at  ERB  sites. 

Much  of  the  organic  carbon  added  to  the  aquifer  is  fermented  to  CH4.  In  coarse-grained 
sediments  near  the  water  table,  CH4  will  degas,  potentially  removing  a  significant  fraction  of  the 
added  electron  equivalents  from  the  saturated  zone.  In  finer  grained  sediments,  gas  bubbles  can 
remain  trapped  within  the  pores  by  capillary  forces.  Once  CH4  production  slows,  the  CH4  in 
these  trapped  bubbles  could  potentially  dissolve  back  into  groundwater  and  be  carried 
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downgradient  with  ambient  groundwater  flow.  CH4  that  is  not  trapped  in  gas  bubbles  or  released 
to  the  vadose  zone  will  migrate  downgradient  with  the  ambient  groundwater  flow.  Aerobic 
methane  oxidation  is  limited  by  a  zone  of  low  O2  and  low  CH4  that  separates  the  anaerobic, 
methane-rich  SWQI  plume  and  the  background  unimpacted,  aerobic  groundwater.  If  AMO 
occurs  using  Fe[III]  as  the  TEA,  then  iron  reduction  will  be  the  main  process  limiting  the  growth 
of  the  dissolved  CH4  plume.  However  if  AMO  rates  are  slow,  dissolved  CH4  could  migrate 
downgradient  with  little  attenuation  beyond  dispersive  mixing.  Once  the  primary  organic 
substrate  (EVO)  is  exhausted,  CH4  production  stops  and  the  CH4  concentrations  will  decline. 

Shortly  after  ED  addition,  SO42'  concentrations  decline  in  the  injection  area  and  a  zone  of 
reduced  SO42"  concentrations  begins  to  migrate  downgradient.  While  sulfide  is  produced  by 
S042"  reduction,  downgradient  migration  is  limited  by  reaction  with  Fe[II]  and  precipitation  as 
low  solubility  sulfide  minerals. 

ED  addition  results  in  release  of  dissolved  Fe2+  with  depletion  of  Fe[III]  in  the  injection  area 
sediments.  However,  downgradient  migration  of  the  dissolved  Fe2+  plume  is  delayed  behind  the 
low  O2  plume  by  sorption  to  the  sediments.  Over  time,  the  zone  of  depleted  sediment  Fe[III] 
grows  is  projected  to  grow  as  dissolved  CH4  reduces  Fe[III]  in  the  downgradient  aquifer.  Once 
the  added  ED  is  degraded,  the  large  reservoir  of  sorbed  Fe[II]  remains  and  is  not  expected  to 
readily  desorb.  Over  time,  dissolved  Fe2+  concentrations  decline  due  to  downgradient  migration 
and  subsequent  sorption  to  the  aquifer  material.  As  clean  aerobic  groundwater  enters  from 
upgradient,  Fe[II]  desorbs  from  the  sediment,  mixes  with  O2,  and  re-precipitates  as  Fe[III]. 

Overall,  Mn  behaves  similarly  to  Fe  with  Mn[IV]  reduced  to  dissolved  Mn2+  in  the  source  area, 
followed  by  downgradient  migration  of  Mn2+.  The  total  amount  of  Mn[II]  produced  is  expected 
to  be  much  lower  than  Fe[II]  due  to  the  smaller  amount  of  Mn[IV]  originally  present  in  many 
aquifer  materials.  The  zone  of  depleted  Mn[IV]  is  projected  to  grow  more  rapidly  than  the 
Fe[III]  depleted  zone  because  Mn  can  be  reduced  under  relatively  more  oxidizing  conditions 
than  Fe.  Similar  to  Fe,  downgradient  migration  of  dissolved  Mn2+  is  limited  by  sorption  to  the 
sediments  and  the  dissolved  Mn2+  dissipates  once  the  ED  is  depleted. 

The  aquifer  downgradient  of  ERB  sites  is  expected  to  be  quickly  depleted  of  O2  following  ED 
addition.  Reoxygenation  of  the  aquifer  is  expected  to  occur  very  slowly,  even  once  the  ED  is 
depleted,  due  to  the  large  reservoir  of  sorbed  Fe[II]  on  the  sediments  that  consumes  O2 
transported  from  upgradient. 


BENEFITS 

This  research  has  improved  our  understanding  of  the  near-  and  long-term  impacts  to  groundwater 
quality  after  implementation  of  in  situ  anaerobic  bioremediation  approaches.  Project  results 
were  used  to  develop  a  general  Conceptual  Model  and  guidance  document  “Extent  and 
Persistence  of  Secondary  Water  Quality  Impacts  after  Enhanced  Reductive  Bioremediation”. 

This  document  provides  information  on  the  microbiological  and  geochemical  processes 
controlling  production  and  attenuation  of  secondary  water  quality  impacts  and  typical 
characteristics  of  secondary  water  quality  impact  plumes  and  can  be  used  by  project  managers, 
consultants,  and  regulators  to  evaluating  the  duration,  extent,  and  magnitude  of  secondary  water 
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quality  impacts  and  to  develop  monitoring  and  management  approaches  to  minimize  these 
impacts. 

Several  results  from  this  study  show  that  natural  attenuation  of  SWQIs  can  be  effective.  The 
first  concerns  the  fate  of  aqueous  Fe[II]  mobilized  during  the  active  phase  of  ERB.  Modeling  of 
the  Bemidji  crude  oil  spill  site  demonstrated  that  91%  of  mobilized  Fe[II]  is  readsorbed  to  the 
sediments  in  the  reducing  zone,  preventing  migration  from  the  treatment  zone  (Ng  et  al.,  2015). 
In  addition,  a  field  study  of  the  fate  of  mobilized  arsenic  at  the  Bemidji  site  shows  that  the 
arsenic  mobilized  by  natural  attenuation  of  crude  oil  is  re-adsorbed  over  the  same  interval  as  the 
Fe[II]  (Cozzarelli  et  al.,  2015).  Production  of  CH4  and  transport  from  the  treatment  zone  can 
also  be  a  concern  during  ERB.  The  Bemidji  modeling  results  show  that  70%  of  the  produced 
CH4  enters  the  gas  phase  and  is  oxidized  before  reaching  the  surface.  This  greatly  limits  the 
migration  of  CH4  in  the  groundwater.  The  recovery  time  of  aerobic  conditions  in  an  aquifer 
following  completion  of  ERB  was  examined  with  modeling  of  the  Cape  Cod  wastewater  study 
site.  The  results  showed  that  recovery  of  aerobic  conditions  can  take  decades,  but  dissolved  Fe2+ 
does  not  appear  to  migrate  long  distances  (Ng  et  al.,  in  prep.). 
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1.0  INTRODUCTION 


1.1  ENHANCED  REDUCTIVE  BIOREMEDIATION  AND  SECONDARY  WATER 
QUALITY  IMPACTS 

Electron  donor  (ED)  addition  can  be  very  effective  in  stimulating  enhanced  reductive 
bioremediation  (ERB)  of  a  broad  range  of  contaminants  including  chloroethenes,  chloroethanes, 
chloromethanes,  chlorinated  cyclic  hydrocarbons,  various  energetics  (e.g.,  perchlorate,  RDX, 
TNT),  and  nitrate  (AFCEE,  NFESC,  ESTCP,  2004).  Hexavalent  chromium  (Cr[VI])  and 
uranium  (U[VI])  can  be  reduced  to  less  mobile  and/or  less  toxic  forms  (e.g.,  Cr[III],  U[IV]) 
(Wielinga  et  al.,  2000;  Yurovsky  et  al.,  2009).  Similarly,  reduction  of  sulfate  to  sulfide 
(sulfidogenesis)  can  be  used  to  precipitate  certain  heavy  metals  (e.g.,  zinc,  cadmium,  and  cobalt) 
(Vanbroekhoven  et  al.,  2009).  However,  ED  addition  can  result  in  the  release  of  a  groundwater 
‘plume’  with  reduced  levels  of  dissolved  oxygen  (O2),  nitrate  (NO3'),  and  sulfate  (SO42'),  and 
elevated  levels  of  dissolved  manganese  (Mn2+),  iron  (Fe2+),  methane  (CH4),  organic  carbon,  and 
naturally  occurring  hazardous  compounds  (e.g.,  arsenic  [As])  (AFCEE,  NFESC,  ESTCP,  2004; 
Hering  et  al.,  2009). 

There  is  growing  concern  about  these  ‘secondary  impacts’  of  ERB.  In  most  cases,  the 
groundwater  downgradient  from  ERB  projects  is  not  toxic,  but  may  have  impaired  taste,  odor, 
and/or  aesthetic  quality,  possibly  making  it  unsuitable  as  a  potable  water  source.  Increasingly, 
regulators  are  requiring  documentation  of  the  expected  or  potential  impacts  prior  to  issuing 
required  permits.  While  secondary  impacts  are  monitored  at  many  ERB  sites,  little  of  this  data 
has  been  compiled  in  a  usable  form  that  can  be  shared  with  regulators.  Fortunately,  there  is  a 
wealth  of  information  on  the  natural  attenuation  of  these  same  parameters  at  sites  where  organic 
rich  materials  have  entered  the  subsurface  (e.g.  wastewater,  landfill,  and  hydrocarbon  plumes). 

For  the  purposes  of  this  report,  we  define  secondary  water  quality  impacts  (SWQIs)  as  changes 
in  groundwater  quality  resulting  from  in  situ  anaerobic  bioremediation  that  have  significant 
adverse  impacts  on  the  beneficial  use  of  the  groundwater  resource.  SWQIs  may  be  temporary  or 
long-lived,  expected  or  unexpected,  have  significant  health  consequences  or  only  result  in  a  mild 
change  in  the  palatability  of  the  water  supply.  SWQIs  do  not  include  contaminants  produced 
from  destruction  of  the  primary  contaminant,  but  are  strictly  limited  to  those  constituents 
produced  by  or  mobilized  due  to  conditions  caused  by  the  in  situ  technology. 

1.2  TECHNICAL  OBJECTIVES 

The  overall  objective  of  this  research  is  to  develop  an  improved  understanding  of  the  near-  and 
long-term  impacts  to  groundwater  quality  after  implementation  of  in  situ  anaerobic 
bioremediation  processes.  Specific  objectives  achieved  during  this  project  include: 

A.  Formulate  a  general  modeling  approach  appropriate  for  simulating  the  production  and 

subsequent  natural  attenuation  of  SWQIs  for  a  wide  range  of  sites.  The  model  approach  and 
validation  was  based  on  data  from  two  intensively  studied  field  sites  —  Bemidji  Crude  Oil 
Spill  and  Cape  Cod  Wastewater  Plume. 
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B.  Assemble  a  database  of  SWQIs  at  ERB  sites  to  determine  the  range  of  SWQI  concentrations 
and  evaluate  trends  over  time  and  distance  downgradient. 

C.  Use  the  validated  modeling  approach  and  results  from  the  database  analysis  to  develop  a 
general  conceptual  model  of  SWQI  formation,  mobilization,  and  attenuation.  This 
conceptual  model  can  be  used  to  identify  sites  where  SWQI  could  be  an  issue  and  develop 
monitoring  and  management  approaches  to  minimize  adverse  impacts. 

1.3  REPORT  ORGANIZATION 

This  technical  report  describes  work  completed  to  date  on  this  project.  Section  2  describes 
model  development  and  evaluation  at  the  Bemidji  and  Cape  Cod  sites.  Section  3  describes 
development  of  the  SWQI  database  and  statistical  analysis  results.  Section  4  presents  the 
indicator  simulation  results.  Section  5  integrates  the  results  from  Sections  2,  3  and  4  to  develop 
a  general  conceptual  model  of  SWQI  production  and  attenuation.  Research  implications  and 
opportunities  are  presented  in  Section  6. 
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2.0  SWQI  MODEL  DEVELOPMENT 


2.1  INTRODUCTION 

One  of  the  goals  for  this  project  was  to  formulate  a  reactive  transport  model  to  simulate  natural 
attenuation  of  the  most  problematic  chemical  species  driving  SWQIs  at  ERB  sites.  These 
species  include  Fe2+,  Mn2+,  sulfide  (S2'),  CH4,  organic  carbon,  and  total  dissolved  solids.  The 
behavior  of  these  and  other  reactive  species  in  groundwater  have  been  intensively  studied  for 
many  years  by  the  USGS  Toxics  Substances  Hydrology  Program  (Toxics  Program).  The 
modeling  goals  of  this  task  were  accomplished  using  datasets  from  two  Toxics  Program  sites:  (a) 
Bemidji  Crude  Oil  Spill,  MN;  and  (b)  Cape  Cod  Wastewater  Plume,  MA.  Together  these  sites 
provided  a  comprehensive  dataset  to  constrain  model  formulations  for  the  suite  of  reactions 
controlling  the  fates  of  the  chemical  species  listed  above. 

Data  from  the  three  sites  together  with  existing  model  results  were  used  to  address  multiple 
issues  concerning  appropriate  representations  of  model  reactions.  First,  the  important  chemical 
species  controlling  the  reaction  behavior  of  the  target  species  were  identified.  Probably  the  most 
important  contribution  of  this  modeling  effort  is  advanced  treatment  of  reactions  between 
aqueous  and  solid  species.  Most  of  the  important  reactions  controlling  Fe2+,  Mn2+  are  reactions 
of  this  type.  During  ERB,  reduced  species  are  often  sequestered  on  solids  in  the  treatment  zone 
and  the  downgradient  aquifer,  such  that  reduced  conditions  persist  after  treatment  has  ended.  A 
goal  of  the  modeling  was  to  provide  a  better  understanding  of  the  time  required  for  recovery  of 
aerobic  conditions  after  ERB  has  ended.  At  Cape  Cod,  the  slow  rate  of  aquifer  recovery  since 
sewage  disposal  ceased  provided  an  extreme  end  member  for  recovery  rate. 

The  sections  below  provide  detailed  descriptions  of  the  model  results  from  each  site. 


2.2  THE  BEMIDJI  CRUDE  OIL  PLUME 
2.2.1  Background  and  Objectives 

Any  source  of  reduced  carbon  emplaced  in  the  subsurface  can  create  SWQIs,  including 
petroleum  hydrocarbons  from  leaks  and  spills  of  crude  oil  and  refined  fuels.  In  1979,  a  pipeline 
rupture  resulted  in  infiltration  of  -2,500  barrels  of  crude  oil  to  a  shallow  glacial  outwash  aquifer 
near  Bemidji,  Minnesota.  The  oil  is  present  in  three  oil  bodies  at  the  water  table  3-8  meters  (m) 
below  the  land  surface.  One  of  the  three  oil  bodies,  known  as  the  “north  pool”,  has  been 
intensively  monitored  by  the  USGS  for  over  30  years.  Site  research  activities  include  detailed 
characterization  of  aquifer  microbiology,  geochemistry,  and  mineralogy  (Essaid  et  al.,  2011). 
Important  SWQIs  at  the  site  result  from  CH4  production  (Amos  et  al.,  2005),  formation  of  Fe2+ 
and  Mn2+  (Baedecker  et  al.,  1993),  and  transport  of  non-target  dissolved  organic  carbon  (DOC) 
compounds  (Thom  and  Aiken,  1998). 

The  objective  of  the  study  was  to  evaluate  the  biogeochemical  processes  that  control  the 
spatiotemporal  evolution  of  both  the  primary  organic  carbon  source  and  SWQIs.  The  focus  at 
this  site  was  to  better  understand  the  mechanisms  leading  to  the  formation  and  natural 
attenuation  of  the  Fe2+,  organic  carbon,  and  CH4  plumes.  Model  simulations  (Ng  et  al.  2015) 
were  used  to  explore  different  conceptual  models  of  CH4  generation  and  outgassing,  formation 
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of  sediment-bound  reduced  iron,  and  pH  buffering  postulated  by  Ng  et  al.  (2014). 

2.2.2  Approach 

The  approach  was  to  create  a  comprehensive  reactive  transport  model  that  accounts  for 
oxidation-reduction  (redox)  reactions  driven  by  organic  carbon  and  equilibrium  reactions  with 
aquifer  minerals  and  groundwater.  The  model  accounts  for  all  features  identified  as  important  in 
previous  Bemidji  models  including  different  classes  of  organic  carbon  in  the  oil  source 
(including  benzene,  ethybenzene,  xylene  (BEX);  toluene;  nonvolatile  dissolved  organic  carbon 
(NVDOC),  and  /2-alkanes),  outgassing  of  CH4  and  carbon  dioxide  (CO2),  multiple  pH  buffering 
mechanisms,  Fe  and  Mn  mineral  phases,  reactions  with  carbonates,  and  representation  of 
sorption  of  Fe,  Mn,  and  the  hydrogen  ion  (H+)  (Ng  et  al.  2015).  The  conceptual  model 
implemented  is  based  on  the  data  compilation  and  mass  balance  analyses  described  in  Ng  et  al. 
(2014)  and  illustrated  in  Figure  2.1. 


Figure  2.1.  Conceptual  model  of  the  Bemidji  north  pool  implemented  with 
the  reactive  transport  model.  Included  are  processes  not  considered  in  previous  Bemidji 
modeling  studies,  such  as  distinct  degradation  pathways  for  diverse  oil  constituents,  direct 
outgassing  in  addition  to  outgassing  from  the  aqueous  phase,  and  Fe  sorption  with  H+  exchange. 

Taken  from  Figure  (1)  in  Ng  et  al.  (2015). 


Data  from  30  years  of  monitoring  were  compiled  to  constrain  the  model  (Figure  2.2).  New 
datasets  that  had  not  been  available  in  the  past  were  also  incorporated,  including  oil  composition, 
sediment  iron  extractions,  and  surface  CO2  efflux.  Data  representing  four  time  periods  were 
used  to  constrain  the  model:  (1)  the  time  of  the  oil  spill  in  1979;  (2)  1987  ±  2  years;  (3)  1993  ±  2 
years;  and  (4)  2008  ±  2  years.  The  data  categories  include  (1)  oil  characterization  and 
quantification;  (2)  concentrations  of  DOC;  (3)  aqueous  concentrations  of  inorganic  groundwater 
plume  constituents;  (4)  solid  iron  concentrations  from  sediment  extractions,  and  (5)  surface 
efflux  of  CO2.  A  detailed  list  of  data  sources  is  provided  by  Ng  et  al.  (2014)  and  the  data  used  in 
the  model  are  available  at  http://mn.water.usgs.gov/projects/bemidji/data/. 

The  model  PHT3D  (Prommer  et  al.,  2003)  was  used  to  simulate  plume  growth  for  30  years  along 
a  2-D  cross-section  of  the  Bemidji  north  pool  transect.  Reactions  were  implemented  using  the 
partial  equilibrium  (PE)  approach  (Jakobsen  and  Postma,  1999)  and  microbial  growth  was 
neglected.  In  the  PE  approach,  organic  carbon  oxidation  is  the  rate-limiting  step  and  reduction 
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of  other  constiutents  is  represented  as  equilibrium  reactions.  For  example,  for  one  mole  of 
benzene  (CeHe)  degrading  at  a  specified  rate,  six  moles  of  C(-l)  and  six  moles  of  H(+l)  are 
newly  introduced  into  the  equilibrium  solution  (oxidation  states  in  parentheses).  C(-l)  then 
automatically  converts  to  the  inorganic  carbon  form  via  the  Terminal  Electron  Accepting  Process 
(TEAP)  with  the  most  favorable  equilibrium  thermodynamic  conditions  in  the  solution.  The 
equilibrium  redox  couples  needed  to  simulate  the  secondary  plumes  at  Bemidji  are  H  '/Fh, 
O2/H2O,  Fe[III]/Fe[II],  Mn[IV]/Mn[II],  and  CO2/CH4.  Other  aqueous  phase  equilibrium 
reactions  considered  include  those  found  in  the  PHT3D  database,  such  as  full  carbonate  and 
hydroxide  chemistry  and  re-oxidation.  Tables  listing  the  oil  components,  mineral  reactions, 
initial  and  boundary  conditions,  and  reaction  parameters  are  provided  in  Ng  et  al.  (2015). 
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Figure  2.2.  Observations  of  the  six  major  redox  species  at  the  three  observational  periods. 

Plots  include  well  data  from  the  available  year  closest  to  the  target  year  (within  ±2  years)  and  are 
adapted  to  our  model  resolution.  Left-hand  stripes  indicate  the  background  and  initial 
concentrations,  determined  from  upgradient  well  observations.  The  boxes  correspond  to  the  oil 
body  reference  rectangle.  The  polygons  outline  the  plume  extent  for  the  corresponding 
component,  which  are  delineated  around  wells  with  measurements  significantly  different  than 
background.  Adapted  from  Figures  (5)  and  (6)  from  Ng  et  al.  [2014]. 
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2.2.3  Major  Findings 

The  model  results  for  the  major  redox  species  are  shown  in  Figure  2.3.  The  model  results 
capture  important  features  of  the  plume.  These  include  the  limited  size  of  the  BEX  plume,  the 
growth  of  the  NVDOC  plume,  the  production  and  subsequent  sorption  of  the  Fe[II]  and  Mn[II] 
plumes,  the  production  and  attenuation  of  dissolved  CH4,  and  finally  the  production  and  growth 
of  the  dissolved  inorganic  carbon  (DIC)  plume.  The  most  important  findings  of  the  modeling 
study  are  summarized  below  and  described  in  detail  by  Ng  et  al.  (2015). 

The  first  major  finding  of  the  modeling  study  concerns  the  fate  of  the  reduced  Fe  produced  by 
dissimilatory  iron-reduction.  The  Bemidji  results  show  that  99%  of  the  total  reduced  Fe  is  in  the 
solid  phase  on  the  surface  of  the  aquifer  sediments.  Immobilization  on  sediments  strongly  limits 
the  SWQI  from  reduced  Fe  at  this  site.  The  model  results  suggest  that  the  immobilization 
mechanism  is  sorption.  The  role  of  sorption  in  sequestration  of  reduced  Fe  appears  much  greater 
than  previously  recognized  at  this  and  other  sites.  The  model  also  showed  that  precipitated 
siderite  serves  as  an  additional  Fe  sink  after  existing  Fe  sorption  sites  are  exhausted.  This 
additional  reduced  Fe  sink  can  facilitate  continued  and  gradual  depletion  of  ferric  oxide  minerals. 
The  form  of  reduced  Fe  on  sediments  may  affect  its  stability,  which  controls  the  likelihood  of 
remobilization  and  future  plume  persistence. 

The  second  major  finding  concerns  the  production  and  fate  of  CH4.  Most  of  the  CH4  is  produced 
by  degradation  of  the  //-alkane  fraction  of  the  oil  (Ng  et  al.,  2014).  The  produced  CH4  and  CO2 
migrate  away  from  the  oil  in  both  the  gas-  and  aqueous-phases  (Amos  et  al.,  2005).  Depleted 
nitrogen  (N2)  concentrations  in  wells  downgradient  from  the  oil  body  indicate  some  outgassing 
of  dissolved  CH4  has  occurred  (Amos  et  al.,  2005),  but  dissolved  N2  is  not  completely  depleted, 
providing  a  constraint  on  the  CH4  cycled  through  the  aqueous  phase.  The  model  results  predict 
that  86%  of  observed  surface  efflux  consists  of  outgassed  carbon  that  could  not  have  cycled 
through  the  aqueous  phase  sampled  in  the  monitoring  wells.  Direct  outgassing  is  the  largest  sink 
for  carbon  from  the  source  zone,  with  an  estimated  70%  of  carbon  lost  by  this  pathway. 
Outgassing  from  the  aqueous  phase  is  mostly  in  the  form  of  CH4,  which  is  consistent  with  its 
much  lower  solubility  compared  to  CO2.  Preferential  outgassing  of  dissolved  CH4  maintains 
low  CH4  to  CO2  aqueous  concentration  ratios,  even  when  the  plume  is  strongly  methanogenic. 

Two  important  attenuation  mechanisms  limit  the  migration  of  the  produced  CH4.  The  gas-phase 
CH4  that  diffuses  upward  through  the  vadose  zone  is  oxidized  to  CO2  before  reaching  land 
surface.  The  primary  natural  attenuation  mechanism  for  aqueous-phase  CH4  in  groundwater  at 
the  Bemidji  site  appears  to  be  anaerobic  methane  oxidation  (AMO)  coupled  to  Fe  reduction 
(Amos  et  al.,  2012).  The  existence  of  AMO  coupled  to  iron  reduction  has  yet  to  be  fully 
established  by  identification  of  the  responsible  microorganisms.  Nevertheless,  field  evidence  at 
Bemidji  (Amos  et  al.,  2012)  and  other  sites  indicate  that  AMO  coupled  to  Fe  reduction  may  play 
an  important  role  in  limiting  the  SWQIs  of  CH4  produced  in  ERB  source  zones  at  some  sites. 
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Figure  2.3.  Model  simulation  results  for  six  major  redox  species  at  three  times 
corresponding  to  the  observational  time  periods  chosen  for  the  data  compilation. 

Gray  boxes  correspond  to  the  oil  body  reference  rectangle.  BEX  and  NVDOC  represent  electron 
donors.  Mn2+,  Fe2+,  and  CH4  are  reduced  aqueous  species  resulting  from  different  anaerobic 
TEAPs.  DIC  is  produced  by  all  TEAPs. 


2.2.4  Important  Implications 

The  Bemidji  model  results  suggest  that  sorption  is  the  dominant  immobilization  mechanism. 
Further  studies  of  this  mechanism  described  in  the  Cape  Cod  section  explore  the  role  of  sorbed 
Fe[II]  as  aerobic  conditions  are  re-established  in  an  aquifer. 

Methane  production  occurs  primarily  in  the  source  zone  where  low  solubility  oil  components 
degrade  without  entering  the  groundwater  plume.  The  CH4  is  produced  in  sufficiently  large 
quantities  that  it  enters  the  gas  phase  and  diffuses  toward  the  surface.  Production  of  CH4  in  the 
source  zone  is  likely  to  occur  when  excess  carbon  is  emplaced  at  ERB  sites.  The  risk  of  upward 
migration  of  CH4  toward  the  land  surface  depends  on  the  nature  and  depth  of  the  emplaced 
carbon. 

The  database  survey  shows  that  dissolved  CH4  plumes  extend  more  than  50  m  downgradient  of 
the  injection  zone  at  some  ERB  sites,  while  at  other  sites  the  CH4  plume  appears  to  be  more 
limited.  The  Bemidji  results  provide  evidence  that  AMO  coupled  to  Fe  reduction  may  be  an 
important  natural  attenuation  mechanism  for  dissolved  CH4  plumes  at  some  ERB  sites. 
Additional  research  is  needed  on  this  process  and  its  reliability  at  other  sites. 

2.3  THE  CAPE  COD  WASTEWATER  PLUME 
2.3.1  Background  and  Objectives 

Almost  two  decades  of  groundwater  chemistry  observations  in  a  wastewater  contaminant  plume 
following  removal  of  the  source  have  led  to  new  insights  into  aquifer  recovery  following  source- 
cessation.  The  contaminant  plume  resulted  from  nearly  60  years  of  land-disposal  of  treated 
wastewater  on  western  Cape  Cod,  Massachusetts,  USA.  The  plume  affects  an  unconfined 
aquifer  with  glacial  outwash  sediments  comprised  of  coarse  sand  and  gravel.  High  sediment 
permeability  combined  with  a  humid  climate  result  in  groundwater  velocities  of  approximately 
0.4  meters  per  day  (m/d)  and  rapid  plume  growth  to  a  length  of  over  7000  meters  in  the  direction 
of  groundwater  flow.  The  aquifer  sediments  support  poor  soils  and  have  very  low  concentrations 
of  organic  carbon.  As  a  result  of  low  loads  of  organic  carbon  to  the  aquifer,  uncontaminated 
groundwater  remains  oxic  with  only  a  small  loss  of  dissolved  oxygen  over  several  kilometer- 
long  transport  distances.  In  contrast,  during  active  wastewater  disposal  the  core  of  the 
wastewater  plume  became  anaerobic  at  the  source.  Monitoring  of  the  plume  within  600  meters 
of  the  abandoned  disposal  beds  since  source-cessation  has  shown  that  unreactive  wastewater 
constituents  (e.g.,  boron  [B])  were  flushed  from  the  study  area  within  10  years  but  the  core  of  the 
plume  has  remained  anoxic  with  a  zone  dominated  by  iron(II)  (Fe[II])  (Figure  4.4).  The 
accumulation  of  Fe[II]  was  the  result  of  anaerobic,  microbial  reduction  of  Fe[III]  oxyhydroxides 
associated  with  sediments,  which  is  driven  by  inputs  of  wastewater-derived  organic  carbon  to  the 
aquifer. 
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Figure  2.4.  Spatial  distribution  of  dissolved  constituents  over  time.  Dissolved  boron  (top), 
oxygen  (middle)  and  Fe(II)  along  a  two-dimensional  cross  section  through  the  wastewater  plume 
near  the  time  of  source  cessation  (left  column),  four  years  after  source  cessation  (center  column), 

and  11  years  after  source  cessation  (right  column). 


The  modeling  objectives  of  plume  evolution  at  the  Cape  Cod  site  are  three-fold.  First,  we  want 
to  develop  a  transferable  modeling  approach  to  simulate  redox  evolution  driven  by  microbial 
metabolic  processes.  Second,  we  want  to  illustrate  the  importance  of  accounting  for  sorption  of 
Fe[II],  how  Fe[II]  sorption  is  coupled  to  sorption  of  other  solutes,  and  how  these  sorption 
reactions  are  coupled  to  microbial  redox  reactions.  Third,  we  want  to  understand  the  processes 
controlling  re-oxygenation  of  the  anaerobic  plume  following  source-cessation. 

2.3.2  Approach 

In  order  to  achieve  the  modeling  objectives,  we  first  developed  a  reactive  transport  model  that 
simulated  the  changes  in  O2,  Fe[H],  pH,  and  other  constituents  observed  over  the  nearly  10  years 
since  source-cessation.  We  then  applied  this  model  to  predict  the  evolution  of  redox  conditions 
in  the  plume  until  the  aquifer  became  fully  oxygenated  within  the  spatial  domain  of  the  model. 
The  microbial  metabolism  during  recovery  is  fueled  by  the  labile  fraction  of  organic  carbon  on 
the  aquifer  sediments.  The  source  of  this  labile  organic  carbon  appears  to  be  mainly  from 
microbial  growth  and  microbial  polymers  formed  during  the  period  of  active  wastewater  disposal 
into  the  aquifer.  We  hypothesize  that  this  labile  organic  carbon,  which  will  be  called  biogenic, 
mixed-stoichiometry  organic  carbon  (BMSOC),  is  a  very  small  fraction  of  the  total  organic 
carbon  on  the  sediment  and  may  not  correlate  with  organic  carbon  quantified  using  typical 
assays.  We  used  the  model  to  examine  the  role  of  BMSOC  and  sediment-bound  Fe[II]  in 
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limiting  the  advancement  of  O2  through  the  plume  and  examined  the  sensitivity  of  the  simulated 
evolution  of  plume  chemistry  to  key  parameters  in  the  model. 

The  model  was  constructed  along  a  two-dimensional  vertical  transect  through  the  aquifer.  A 
two-dimensional,  steady,  groundwater  flow  field  was  set  up  to  resemble  the  average  flow  field 
observed  in  the  study  area  and  to  match  the  slope  along  the  direction  of  flow  of  the  vertical 
boundary  between  the  anoxic  plume  and  oxic,  uncontaminated  zone  above  it.  The  observed 
biogeochemical  zones  in  the  plume  prior  to  source  cessation  were  established  within  the  spatial 
domain  of  the  model  and  include:  oxic,  anoxic  with  nitrate,  and  anoxic  with  Fe[II]  and 
manganese  (Mn[III/II])  (Figure  2.5).  A  set  of  chemical  reactions  was  included  to  describe 
competitive  sorption  among  the  major  cations  (Na+,  K+,  Mg2+,  Ca2+),  Fe[II]  and  Mn[III/II],  and 
H+.  Parameters  in  the  sorption  model  were  calibrated  using  laboratory  experiments  with  site- 
specific  sediments  (Kohler  and  Kent,  in  prep.).  The  initial  spatial  distribution  of  BMSOC  and 
the  rates  of  aerobic  respiration  and  anaerobic  respiration  (coupled  to  NO3'  or  Fe[III]  reduction) 
were  calibrated  using  1)  the  observed  slow  advance  of  the  front  of  O2  beneath  the  abandoned 
wastewater  disposal  beds;  2)  the  observed  evolution  of  redox  conditions  within  the  plume 
downgradient  of  the  source  following  source  cessation;  and  3)  the  evolution  of  pH  in  the 
anaerobic  plume  following  source  cessation. 
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Figure  2.5.  Groundwater  flow  field  (arrows)  and  geochemical  zones  in  groundwater  model. 

Biogeochemical  zones:  I  (oxic,  uncontaminated);  II  (oxic,  contaminated,  directly  under  disposal 
beds;  III  (anoxic  with  nitrate,  under  disposal  beds  and  immediately  downgradient);  IV  (anoxic 
with  nitrate);  and  anoxic  with  Fe[II]  (Va)  or  Mn  (Vb)  along  the  two-dimensional  vertical  cross 

section  model  domain. 


2.3.3  Major  Findings 

Simulations  show  that  the  O2  front  advances  much  more  slowly  than  the  groundwater  velocity 
(Ng  et  al.,  in  prep.),  as  was  observed  in  the  field  (Figure  2.4).  Unreactive  constituents  were 
transported  out  of  the  model  domain  in  a  few  years  but  the  model  predicts  that  re-oxygenation 
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will  require  a  few  decades.  Simulation  results  show  that  both  BMSOC  and  sorbed  Fe[II]  control 
the  time  required  for  re-oxygenation. 

Simulation  results  show  that  BMSOC  is  an  important  control  on  the  rate  of  re-oxygenation  of  the 
wastewater  plume,  especially  in  the  regions  under  the  abandoned  wastewater  disposal  beds  and 
between  the  abandoned  disposal  beds  and  the  Fe[II]  zone  (zone  III  in  Figure  2.5)  (Ng  et  al.,  in 
prep.).  The  distribution  of  BMSOC  in  the  model  is  shown  in  Figure  2.6.  The  maximum 
BMSOC  sediment  concentration,  which  occurs  under  the  abandoned  wastewater  disposal  beds,  is 
0.6  micromoles  C  per  gram  of  sediment  (pmol/g)  (Ng  et  al.,  in  prep.).  For  comparison,  total  and 
water-extractable  organic  carbon  concentrations  observed  on  sediments  under  the  disposal  beds 
were  20-50  and  9-30  pmol/g,  respectively  (Smith  et  al.,  2013).  Previous  calculations  showed 
that  a  value  of  0.7  pmol/g  for  metabolizable  organic  carbon  on  the  sediment  could  account  for 
the  observed  rate  of  travel  of  Oi  under  the  abandoned  disposal  beds  (Smith  et  al.,  2013).  Thus, 
the  concentration  of  sediment-bound  organic  carbon  that  drives  microbial  respiration  in  the 
wastewater  plume  is  very  low  and  does  not  correspond  to  any  form  of  sediment-bound  organic 
carbon  thus  far  measured.  The  model-calibrated  first-order  rate  of  oxygen  consumption  coupled 
to  BMSOC  degradation  (Ng  et  al.,  in  prep.)  is  also  considerably  lower  than  oxygen  consumption 
rates  measured  in  batch  incubations  reported  in  Smith  et  al.  (2013). 
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Figure  2.6.  Initial  spatial  distribution  of  concentrations  of  reactive  organic  carbon  (top), 
dissolved  Fe[II]  (middle),  and  sorbed  Fe[II]  (bottom)  used  in  simulations. 


In  addition  to  BMSOC,  sorbed  Fe[II]  plays  an  important  role  in  delaying  re-oxygenation  of  the 
anoxic  zone  of  the  wastewater  plume  (Ng  et  al.,  in  prep.)  (Figure  2.7).  Sorbed  Fe[II]  owes  its 
existence  to  wastewater  derived  organic  carbon,  which  drove  dissimilatory  Fe[HI]  reduction  in 
the  anoxic  core  of  the  plume.  Observed  growth  of  the  Fe[II]  zone  after  cessation  of  wastewater 
disposal  (Figure  2.4),  driven  by  anaerobic  microbial  metabolism,  was  used  to  calibrate  the  initial 
concentration  of  BMSOC  in  the  region  downgradient  from  the  abandoned  disposal  beds  and  to 
constrain  the  rate  constant  for  anaerobic  biodegradation.  The  relationship  between  dissolved  and 
sorbed  Fe[II]  concentrations  is  described  by  the  sorption  model  and  controlled  by  pH  and  the 
major  cation  concentrations.  The  model-predicted  concentrations  of  sorbed  Fe[II]  in  the  Fe[II] 
zone  at  the  time  of  source  cessation  is  in  good  agreement  with  extractions  of  Fe[II]  from 
sediments  collected  from  that  zone.  Sorbed  Fe[II]  concentrations  exceed  dissolved 
concentrations  by  a  factor  of  20  or  more,  depending  on  local  chemical  conditions. 
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Figure  2.7.  Simulated  spatial  distributions  of  dissolved  oxygen  4years  after  source- 
cessation  (left-hand  column),  13  years  after  source-cessation  (center  column),  and  26  years 
after  source-cessation  (right-hand  column)  based  on  three  different  scenarios.  Top  row: 
oxygen  removal  by  oxidation  of  dissolved  and  sorbed  Fe[II]  but  no  removal  coupled  to  microbial 
oxidation  of  organic  carbon.  Middle  row:  oxygen  removal  by  microbial  oxidation  of  organic 
carbon  and  dissolved  Fe[II]  but  no  sorbed  Fe[II].  Bottom  row:  oxygen  removal  by  microbial 
oxidation  of  organic  carbon  and  oxidation  of  dissolved  and  sorbed  Fe[II]. 


The  role  of  aerobic  respiration  coupled  to  BMSOC  and  oxidation  of  sorbed  Fe[II]  in  controlling 
the  rate  of  re-oxygenation  of  the  wastewater  plume  is  compared  in  Figure  2.7.  In  the  absence  of 
aerobic  respiration  coupled  to  BMSOC,  the  dissolved  oxygen  front  advances  rapidly  to  the 
region  with  elevated  sorbed  Fe[II]  (Figure  2.7,  top  row,  left  panel).  Oxidation  of  the  sorbed 
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Fe[II]  is  complete  sometime  prior  to  13  years  after  source  cessation,  at  which  point  the  former 
wastewater  plume  is  predicted  to  be  fully  oxic  (Figure  2.7,  top  row,  middle  panel).  In  the 
absence  of  sorbed  Fe[II]  but  presence  of  dissolved  Fe[II]  at  concentrations  observed  at  the  time 
of  source-cessation  (Figure  2.4),  the  dissolved  oxygen  front  moves  slowly  through  the  aquifer  as 
observed  (Figure  2.7,  center  row).  However,  simulated  pH  values  in  the  downgradient  region 
with  elevated  dissolved  Fe[II]  concentrations  increase  to  the  point  where  the  Fe[II]  carbonate 
mineral  siderite  precipitates  (Ng  et  al.,  in  prep.).  Precipitation  of  siderite  in  the  model  is 
inconsistent  with  field  observations  showing  that  siderite  has  remained  undersaturated  by  a  factor 
of  10  or  more  in  the  wastewater  plume.  The  O2  front  moves  slowly  once  it  encounters  the  region 
with  siderite  (Figure  2.7,  center  row,  right-hand  panel).  In  the  full  simulation  (Figure  2.7, 
bottom  row),  the  O2  front  moves  slowly  directly  under  the  abandoned  wastewater  disposal  beds, 
more  rapidly  in  the  region  between  the  disposal  beds  and  the  region  with  elevated  dissolved  and 
sorbed  Fe[II],  then  slowly  through  the  region  with  sorbed  Fe[II]. 

H+  is  a  reactant  or  product  in  both  redox  (e.g.,  Abrams  et  al.,  1996;  Hunter  et  al.,  1998;  Bethke  et 
al.,  2011)  and  sorption  reactions  (Davis  et  al.,  1998;  Kent  et  al.,  2007).  Therefore,  spatial  and 
temporal  distributions  of  pH  can  provide  important  constraints  for  models  incorporating 
biodegradation.  This  is  especially  true  where  dissimilatory  Fe[III]  reduction  is  important 
because  it  drives  pH  up  (e.g.,  Abrams  et  al.,  1996).  In  our  model,  if  the  rate  of  anaerobic 
biodegradation  is  too  high,  simulated  pH  values  greatly  exceed  those  observed  in  the  field  (Ng  et 
al.,  in  prep.).  Conversely,  if  the  rate  of  anaerobic  biodegradation  is  too  low,  dissolved  Fe[II] 
concentrations  are  unrealistically  low.  Concentrations  of  Fe[II]  and  pH  values  also  constrain  the 
concentration  of  BMSOC  downgradient  of  the  abandoned  wastewater  disposal  beds.  Sorption 
reactions  also  influence  Fe[II]  concentrations  and  pH  values  but  parameters  describing  these 
reactions  are  constrained  by  laboratory  experimental  data  (Kohler  and  Kent,  in  prep.). 
Perturbations  in  sorption  parameters  required  to  mitigate  excessive  pH  values  or  insufficient 
Fe[II]  concentrations  in  the  model  fall  outside  of  the  range  consistent  with  the  results  of 
laboratory  experiments. 

2.3.4  Important  Implications 

Addition  of  the  organic  carbon  source  in  ERB  will  promote  microbial  growth,  which  will 
subsequently  result  in  the  accumulation  of  biogenic  organic  carbon  on  sediments  near  the  source. 
This  biogenic  organic  carbon  may  drive  microbial  metabolism  long  after  depletion  of  the 
injected  carbon  source,  contributing  to  the  persistence  of  anaerobic  conditions  in  the  aquifer  near 
and  downgradient  from  the  source  after  treatment  ends. 

Addition  of  organic  carbon  will  stimulate  production  of  Fe[II]  through  dissimilatory  Fe[III] 
reduction.  The  extent  to  which  sorbed  Fe[II]  accumulates  on  sediments  will  depend  on  local 
chemical  conditions,  namely  pH  and  major  cation  concentrations.  Elevated  pH  values  and  low 
dissolved  salt  concentrations  favor  increased  sorption  and,  therefore  higher  concentrations  of 
sorbed  Fe[II].  Fe [II] -containing  secondary  phases,  such  as  siderite,  green  rust,  or  magnetite, 
could  form  as  well,  but  nucleation  and  growth  of  these  phases  can  be  limited.  Sorbed  Fe[II]  will 
be  an  important  component  of  the  Fe[II]  pool  that  retards  the  re-oxygenation  of  aquifers 
following  ERB  treatments. 

Iron(II)  sorption  is  influenced  by  pH,  which,  in  turn,  is  controlled  by  redox  and  sorption 
reactions.  Reduction  of  Fe[III]  oxyhydroxides  to  Fe[II]  has  an  especially  large  impact  on  pH 
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because  three  moles  of  hydroxide  (OH  )  are  released  per  mole  of  Fe[III]  reduced.  Both  products 
of  Fe[III]  reduction  (OH  and  Fe[II])  participate  in  sorption  reactions  that  also  involve  major 
cations  and  other  groundwater  constituents.  Accounting  for  the  evolution  of  pH  in  predictive 
simulations  requires  proper  description  of  both  redox  and  sorption  reactions.  Because  of  this, 
observations  of  pH  can  provide  insight  into  important  process  occurring  in  the  sub-surface  and 
can  help  constrain  models  for  predicting  secondary  water  quality  impacts  of  ERB. 
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3.0  SWQI  DATABASE 


3.1  INTRODUCTION 

While  ED  addition  will  drive  the  biotic  and  abiotic  destruction  and/or  immobilization  of  a  broad 
range  of  groundwater  contmainants,  it  often  results  in  the  production  of  SWQIs  at  ERB  sites. 
These  SWQIs  could  potentially  have  adverse  impacts  on  the  beneficial  use  of  the  groundwater 
resource.  SWQIs  may  be  temporary  or  long-lived,  expected  or  unexpected,  have  significant 
health  consequences  or  only  result  in  a  mild  change  in  the  palatability  of  the  water  supply.  For 
this  article,  SWQIs  do  not  include  degradation  products  produced  from  the  primary  contaminant, 
but  are  limited  to  those  constituents  produced  by  or  mobilized  due  to  conditions  caused  by  ERB. 

The  most  common  approach  for  stimulating  ERB  involves  addition  of  a  fermentable  organic 
substrate  to  the  aquifer  to  serve  as  an  electron  donor  and  carbon  source  for  microbial  growth. 

The  choice  of  substrate  and  the  method  of  injection  depend  on  the  contaminant  type  and 
distribution  in  the  aquifer,  hydrogeology,  and  remediation  objectives.  Organic  substrate  can  be 
added  using  conventional  well  installations,  by  direct-push  technology,  or  by  excavation  and 
backfill.  Injection  points  may  be  arranged  in  a  grid  to  treat  source  areas  or  in  a  line 
perpendicular  to  groundwater  flow  to  form  a  PRB.  Slow-release  products  composed  of  edible 
oils  or  solid  substrates  tend  to  stay  in  place  for  an  extended  treatment  period.  Soluble  substrates 
or  soluble  fermentation  products  of  slow-release  substrates  can  potentially  migrate  into  and 
disperse  within  heterogeneous  lithologies  via  advection  and  diffusion,  providing  broader  but 
shorter-lived  treatment  zones.  The  added  organic  substrates  are  first  fermented  to  hydrogen  (H2) 
and  volatile  fatty  acids  (VFAs).  The  VFAs,  including  acetate,  lactate,  propionate,  and  butyrate, 
provide  carbon  and  energy  for  bacterial  metabolism. 

At  a  few  sites,  concerns  have  been  raised  regarding  the  formation  and  mobilization  of  SWQIs 
following  ERB.  USEPA  (2009)  has  established  secondary  maximum  contaminant  levels 
(MCLs)  for  a  number  of  SWQIs  including  Fe,  Mn,  and  total  dissolved  solids.  Hering  et  al. 
(2009)  and  He  et  al.  (2010)  showed  that  arsenic  could  be  temporarily  mobilized  during  ERB  of 
chlorinated  solvents  at  Fort  Devens,  Massachusetts.  At  the  Grants  Chlorinated  Solvent  Plume 
Superfund  Site,  a  vapor  intrusion  program  was  implemented  to  determine  if  there  was  potential 
for  methane  and  hydrogen  sulfide  accumulation  (USEPA  2013). 

To  gain  a  better  understanding  of  the  formation  and  extent  of  SWQIs  at  ERB  sites,  monitoring 
data  from  47  active  and  former  ERB  sites  was  compiled  and  analyzed  to  determine  the  range  of 
SWQIs  and  evaluate  trends  over  time  and  distance  downgradient.  This  information  was  then 
used  to  develop  a  general  conceptual  model  of  SWQI  production  in  treatment  zones  and 
attenuation  during  downgradient  transport. 
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3.2  DATABASE  DEVELOPMENT 

SWQI  data  were  collected  from  47  ERB  sites  across  the  United  States.  These  sites  are 
illustrative  of  the  wide  range  of  conditions  found  at  ERB  field  sites,  including  different 
contaminants  of  concern  (COCs),  injection  strategy,  organic  substrates,  and  geographic  location. 
Data  was  collected  from  regulatory  monitoring  reports,  DOD  (i.e.,  AFCEE,  ESTCP,  SERDP) 
project  reports,  journal  articles,  and  personal  correspondence  with  site  project  managers.  A  list 
of  the  sites  included  in  the  database,  including  site  name,  location,  injection  type,  organic 
substrate,  data  source,  and  length  of  data  record,  is  provided  in  Table  A-l  in  Appendix  A. 

Once  site  reports  containing  pre-  and  post-injection  SWQI  data  were  collected,  the  data  were 
tabulated  for  statistical  analysis.  Data  collected  from  each  site  included:  concentrations  of  COCs 
and  daughter  products,  SWQIs,  and  terminal  electron  acceptors  (TEAs);  field  parameter  values 
(e.g.,  pH,  specific  conductance,  turbidity);  injection  data  (injection  strategy,  amount  and  type  of 
ED,  addition  of  other  amendments);  and  hydrogeologic  information  (groundwater  velocity, 
hydraulic  conductivity,  lithology).  Concentrations  reported  as  below  detection  were  entered  into 
the  database  at  half  the  detection  limit,  with  a  note  stating  the  parameter  was  not  detected. 
Summary  statistics  (mean,  median,  standard  deviation,  etc.)  and  cumulative  frequency 
distribution  plots  were  generated  for  upgradient,  treatment  zone,  and  various  downgradient 
locations. 

3.3  REMEDIATION  SYSTEM  CHARACTERISTICS 

The  ERB  projects  included  in  our  survey  used  a  variety  of  different  organic  substrates  including 
soluble  substrates  (14  sites),  low  solubility  liquid  substrates  (21  sites),  solid  substrates  (2  sites), 
and  mixed  or  multiple  substrates  (10  sites).  The  injected  substrates  were  introduced  into  the 
subsurface  as  source  area  point  injections  (32  times),  PRBs  (16  times),  recirculation  systems  (7 
times),  biowalls  (mulch  or  mulch  and  soybean  oil)  (4  times),  or  horizontal  wells  (1  time).  The 
sum  of  these  is  greater  than  47  because  some  sites  received  multiple  injections  as  different 
injection  types. 

In  addition  to  the  primary  substrates,  some  ERB  sites  were  also  amended  with  pH  buffering 
reagents  (e.g.,  bicarbonate,  carbonate)  (10  sites);  bioaugmentation  cultures  (e.g.,  KB-1®,  SDC-9) 
(9  sites);  in  situ  chemical  reduction  reagents  (e.g.,  ferrous  sulfate,  gypsum)  (5  sites);  and 
nutrients  to  aid  biological  growth  (e.g.,  diammonium  phosphate,  yeast  extract,  vitamin  B 12)  (5 
sites).  Available  monitoring  data  ranged  from  4  months  to  1 1  years  post-injection.  When 
additional  iron  (5  sites)  or  sulfate  (2  sites)  was  added  as  part  of  the  remediation  process,  data  on 
iron  and/or  sulfate  concentrations  from  these  sites  were  not  included  in  the  general  statistical 
analysis. 

Table  3.1  shows  the  range  of  remediation  site  characteristics  including  number  of  monitoring 
wells,  monitoring  duration,  estimated  groundwater  velocity  at  the  site,  vertical  injection  interval, 
width  of  barrier  or  area  treatment  perpendicular  to  flow,  length  of  treatment  zone  parallel  to 
flow,  and  amount  of  carbon  added.  Estimated  groundwater  velocity  varied  greatly  from  2  to  350 
meters  per  year  (m/yr).  Many  of  the  remediation  systems  were  relatively  small  and  served  as 
pilot  systems  to  evaluate  system  performance  before  moving  to  full-scale  remediation.  The 
median  length  of  the  injection  zone  parallel  to  flow  was  small  since  many  of  the  remediation 
systems  were  designed  to  operate  as  PRBs  to  limit  contaminant  migration. 
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Table  3.1.  Treatment  System  Characteristics  in  Survey  Population. 


Parameter 

Units 

n 

Percenti] 

Les 

Min. 

10th 

25th 

50th 

75th 

90th 

Max. 

Monitoring  wells  per  site 

wells 

47 

2 

6 

9 

14 

24 

40 

136 

Monitoring  Duration 

yr 

47 

4 

7 

12 

22 

50 

77 

132 

GW  Velocity 

m/yr 

36 

0.8 

3 

5 

35 

94 

250 

350 

Travel  time  to  most 
downgradient  well 

yr 

36 

0.03 

0.06 

0.2 

0.6 

2.7 

18 

48 

Vertical  Treatment  Interval 

m 

33 

2 

3 

3 

6 

6 

9 

20 

Monitoring  Period  Post- 
Injection 

yr 

47 

0.3 

0.6 

1.0 

1.8 

4.1 

6.4 

11 

Monitoring  Period/  Travel 
time 

yr/yr 

36 

0.02 

0.1 

0.5 

2.5 

10 

17 

43 

Treatment  Width 
Perpendicular  to  GW  Flow 

m 

33 

1 

4 

9 

14 

30 

67 

390 

Treatment  Length  Parrallel 
to  GW  Flow 

m 

33 

0.5 

1 

3 

3 

9 

15 

76 

Carbon  Injected 

Tonne 

33 

0.03 

0.12 

0.33 

1.1 

2.7 

11 

70 

Injection  Barrier  Carbon 

kg/m2 

13 

2 

5 

7 

19 

30 

61 

129 

Mulch  Biowall  Carbon 

kg/m3 

2 

19 

19 

20 

21 

21 

22 

22 

Area  Treatment  Carbon 

kg/m3 

18 

0.2 

0.4 

1 

3 

8 

33 

72 

3.4  SWQI  PARAMETER  CUMULATIVE  FREQUENCY  DISTRIBUTIONS 

Table  A-2  in  Appendix  A  provides  summary  statistics  for  post-injection  monitoring  data  for  all 
upgradient,  injection  area,  and  downgradient  wells  at  the  47  sites  contained  in  the  database. 
Downgradient  wells  include  monitoring  points  outside  of  the  target  treatment  zone  and  up  to  760 
m  downgradient.  Figure  3.1  shows  cummulative  frequency  distributions  for  most  of  the 
tablulated  SWQI  parameters  in  monitoring  wells  located  upgradient,  within  the  injection  area, 
and  at  varying  distances  downgradient.  USEPA  MCLs  are  shown  for  reference  where  applicable 
(USEPA  2009).  For  parameters  that  generally  increase  following  ERB  (total  organic  carbon 
[TOC],  Fe,  Mn,  As,  S2",  CEU),  maximum  values  are  shown,  while  minimum  values  are  shown 
for  parameters  that  generally  decrease  (oxidation-reduction  potential  [ORP],  SO42").  Maximum 
or  minimum  values  are  reported  because  average  values  may  be  biased  by  measurements 
collected  before  the  organic  substrate  reached  the  well  or  after  it  was  depleted.  For  parameters 
that  may  increase  or  decrease  (pH),  median  values  are  shown.  Figure  3.2  shows  cumulative 
frequency  distributions  for  the  change  in  measured  values  for  these  same  parameters.  Values 
plotted  are  delta  (A)  concentration  maximum  or  minimum  value  (C)  minus  pre-injection  value 
(Co)  or  C/Co,  depending  on  which  approach  best  illustrated  the  results. 
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Figure  3.1.  Cumulative  frequency  distributions  for  SWQI  concentrations  at  47  ERB  sites. 

Figures  show  measured  values  of  upgradient,  injection  area,  and  downgradient  monitoring  wells 
for:  a)  maximum  total  organic  carbon  (TOC);  b)  median  pH  for  non -buffered  sites;  c)  minimum 
dissolved  sulfate  (SO42');  d)  median  pH  buffered  sites;  e)  maximum  dissolved  sulfide  (S2~);  f) 
maximum  dissolved  iron  (Fe);  g)  maximum  dissolved  methane  (CH4);  h)  maximum  dissolved 
manganese  (Mn);  and  i)  maximum  dissolved  arsenic  (As).  MCL  is  shown  by  black  vertical  line. 
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Figure  3.2.  Cumulative  frequency  distributions  for  change  in  SWQI  concentrations  at  47 
ERB  sites.  Figures  show  change  in  measured  values  in  upgradient,  injection  area,  and 
downgradient  monitoring  wells:  a)  total  organic  carbon  (TOC);  b)  pH  for  non-buff ered  sites;  c) 
dissolved  sulfate  (SO42');  d)  pH  for  buffered  sites;  e)  dissolved  sulfide  (S2");  f)  dissolved  iron 
(Fe);  g)  dissolved  methane  (CH4);  h)  dissolved  manganese  (Mn);  and  i)  dissolved  arsenic  (As). 
A  indicates  maximum  observed  value  minus  pre-injection  value  at  that  well.  C/C0  indicates 
maximum  (or  minimum)  observed  value  divided  by  pre-injection  value  at  that  well.  For  pH, 
change  is  difference  between  median  post-injection  and  pre-injection  value  at  the  well. 
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Maximum  observed  TOC  concentrations  and  changes  in  concentrations  are  shown  in  Figures 
3.1a  and  3.2a.  As  expected  following  organic  substrate  addition,  injection  area  concentrations 
were  substantially  elevated  compared  to  upgradient  concentrations.  For  instance,  the  median 
injection  area  TOC  concentration  is  nearly  two  orders  of  magnitude  higher  than  the  median 
upgradient  TOC  concentration.  However,  TOC  concentrations  decline  rapidly  with  distance, 
dropping  almost  an  order  of  magnitude  in  the  0-10  m  wells.  At  greater  than  10  m  downgradient, 
the  median  TOC  was  2-  2.5  times  the  pre-injection  value. 

Figure  3.1b  and  3.1d  shows  median  post-injection  pH  values  for  wells  at  sites  with  and  without 
buffer  addition.  The  ApH  values  shown  in  Figures  3.2b  and  3.2d  are  the  difference  between  the 
median  pH  and  the  value  prior  to  injection.  At  non-buffered  sites,  the  background  pH  is  slightly 
below  neutral  (Figure  3.1b).  The  largest  pH  declines  occur  in  wells  closest  to  the  injection  area 
and  lessen  with  distance  downgradient  (Figure  3.2b).  The  range  of  ApH  values  for  wells  0-10, 
10-25,  25-50,  and  >50  m  downgradient  are  fairly  similar,  indicating  large  pH  changes  are 
typically  restricted  to  the  injection  area.  At  buffered  sites,  the  upgradient/background  pH  values 
are  low,  which  is  likely  why  buffer  is  added  (Figure  3.1d).  At  buffered  sites,  the  largest  pH 
increases  (positive  ApH)  occurs  in  wells  close  to  the  injection  area,  with  smaller  increases  farther 
downgradient  (Figure  3.2d). 

The  concentration  of  aqueous  electron  acceptors  O2  and  NO3"  were  relatively  low  in  upgradient 
wells  (median  =  0.3  mg/L  O2  and  0.1  mg/L  NO3-N)  and  consequently  the  decline  in  these 
parameters  was  small  (Table  A-2).  There  was  little  or  no  rebound  in  these  parameters 
downgradient  of  the  injection  zone,  presumably  due  to  the  low  background  concentrations. 

Minimum  upgradient  (background)  SO42  concentrations  varied  widely  from  <1  mg/L  to  near  the 
gypsum  solubility  limit  (Figure  3.1c).  In  the  injection  area,  SO42  concentrations  were  reduced, 
typically  by  about  1  order  of  magnitude.  However,  SO42'  was  not  completely  depleted  at  most 
sites.  Further  downgradient,  the  change  in  SO42  concentration  was  smaller,  indicating 
substantial  attenuation  of  the  low  SCL2"  plume  (Figure  3.2c).  SO42  reduction  did  result  in 
detectable  increases  in  S2"  concentrations  in  the  injection  area  and  downgradient  wells  (Figure 
3.1e).  However,  S2‘  was  only  monitored  in  about  340  wells  or  35%  of  the  database  wells,  so  care 
should  be  taken  in  extrapolating  these  results  to  other  sites.  In  those  wells  where  sulfide  was 
monitored,  the  increase  in  sulfide  was  typically  less  than  1  mg/L  (Figure  3.2e),  although  sulfide 
increased  by  more  than  10  mg/L  in  a  few  injection  area  wells. 

Maximum  dissolved  Fe,  Mn,  and  arsenic  concentrations  observed  in  upgradient,  injection  area, 
and  downgradient  wells  are  shown  in  Figures  3.1f,  3.1h,  and  3.1i.  The  maximum  changes  in  Fe, 
Mn,  and  arsenic  in  individual  wells  are  shown  in  Figures  3.2f,  3.2h,  and  3.2i.  The  USEPA 
MCLs  for  these  parameters  are  0.3  mg/L  Fe,  0.05  mg/L  Mn  and  0.01  mg/L  arsenic.  Substrate 
addition  during  ERB  results  in  large  increases  in  Fe  concentrations  in  injection  area  wells,  with 
smaller  increases  in  downgradient  wells  (Figure  3.2f).  At  >50  m  downgradient,  the  median 
increase  in  dissolved  Fe  is  less  than  1  mg/L.  This  indicates  that  iron  is  attenuating  during 
downgradient  migration. 

Manganese  concentrations  exceeded  the  secondary  MCL  at  80%  (48  of  61)  of  upgradient  wells, 
94%  (118  of  125)  injection  area  wells,  and  76%  (152  of  201)  of  downgradient  wells  (Figure 
3.1h).  Manganese  concentrations  in  injection  area  wells  were  significantly  elevated  compared  to 
upgradient  concentrations.  At  the  10th  through  90th  percentiles,  the  injection  area  concentration 
was  at  least  one  order  of  magnitude  higher  than  the  corresponding  upgradient  concentration.  The 
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range  of  manganese  concentrations  generally  lessens  with  increasing  distance  downgradient. 
However,  the  highest  downgradient  manganese  concentrations  were  fairly  similar  to  the  highest 
injection  area  manganese  concentrations,  indicating  significant  downgradient  transport  of 
dissolved  manganese  at  some  sites. 

Arsenic  was  only  monitored  in  177  wells  or  19%  of  the  database  wells,  so  care  should  be  taken 
in  extrapolating  these  results  to  other  sites.  Arsenic  concentrations  at  injection  area  wells  were 
elevated  compared  to  upgradient  concentrations  (Figure  3-li)  indicating  ED  addition  often 
results  in  an  increase  in  injection  area  arsenic  concentrations.  The  concentration  ranges  for  the 
downgradient  divisions  are  fairly  similar  to  each  other,  though  slightly  decreasing  with  distance. 
For  instance,  the  50th  percentile  concentration  is  0.014  mg/L  for  wells  0-10  m  downgradient; 

0.01 1  mg/L  for  wells  10-25  m  downgradient;  0.0056  mg/L  for  wells  25-50  m  downgradient;  and 
0.0056  mg/L  for  wells  greater  than  50  m  downgradient.  While  the  median  arsenic  concentrations 
for  wells  25-50  m  and  greater  than  50  m  were  lower  than  the  MCL,  a  greater  percentage  of  these 
wells  exceeded  the  MCL  (33%  and  40%,  respectively)  than  for  upgradient  wells  (10%).  These 
data  indicate  that  arsenic  mobilization  could  be  a  problem  at  distances  greater  than  50  m 
downgradient.  However,  it  should  be  cautioned  that  there  are  very  few  monitoring  points  greater 
than  25  m  downgradient  (16  in  total).  In  addition,  all  wells  greater  than  50  m  downgradient  are 
from  a  single  site  (Ft  Devens,  MA).  Data  from  other  sites  would  help  verify  if  arsenic 
mobilization  is  an  important  concern.  It  should  be  noted  that  there  is  likely  some  selection  bias 
in  sites  with  arsenic  data;  that  is,  arsenic  is  most  often  monitored  at  sites  where  naturally 
occurring  arsenic  is  known  or  expected  to  be  present. 

Figure  3.3  shows  dissolved  arsenic  and  dissolved  iron  concentrations  in  the  same  samples 
collected  from  injection  area  and  downgradient  monitoring  wells,  when  concentrations  of  both 
Fe  and  arsenic  were  above  the  reporting  limit.  Increasing  dissolved  arsenic  concentrations  were 
highly  correlated  with  increasing  dissolved  Fe  concentrations  (p<0.01),  consistent  with  arsenic 
release  during  reduction  of  Fe  oxides  and/or  arsenic  removal  during  sorption/co-precipitation 
with  Fe.  At  sites  sampled  for  both  arsenic  and  iron,  very  high  dissolved  Fe  concentrations  (i.e., 
greater  than  200  mg/L)  resulted  in  arsenic  consistently  above  the  MCL,  while  very  low  dissolved 
Fe  concentrations  (i.e.,  less  than  0.5  mg/L)  resulted  in  arsenic  consistently  below  the  MCL.  As 
mentioned  above,  care  should  be  taken  when  extrapolating  this  relationship  since  arsenic  was  not 
monitored  at  all  sites,  likely  due  to  selection  bias. 


22 


1 


Iron  (mg/L) 


Iron  (mg/L) 


Figure  3.3.  Arsenic  versus  iron  concentrations  for  injection  area  and  downgradient 

monitoring  wells. 


Maximum  CH4  concentrations  at  each  well  are  shown  in  Figure  3.1g.  The  maximum  change  in 
CH4  concentrations  in  individual  wells  are  shown  in  Figure  3.2g.  Upgradient  CH4 
concentrations  were  substantially  lower  than  injection  area  and  downgradient  concentrations  at 
most  percentiles  due  to  methane  production  following  organic  carbon  addition.  The  maximum 
increase  in  CH4  concentration  declines  with  distance,  indicating  some  attenuation  with  distance. 
However,  there  were  a  significant  number  of  wells  at  greater  than  50  m  downgradient  with  high 
CH4  levels  (>10  mg/L),  indicating  the  potential  for  downgradient  migration  of  dissolved  CH4. 

3.5  POTABLE  WATER  IMPACTS 

While  SWQIs  were  detected  in  downgradient  wells,  the  large  majority  of  these  wells  were  within 
the  primary  contaminant  plume.  At  chlorinated  ethene  sites,  only  31  of  380  downgradient  wells 
(8%)  did  not  contain  at  least  one  chlorinated  ethene  (PCE,  TCE,  cis-  1,2-DCE,  trans-  1,2-DCE, 

1,1 -DCE  ,  or  vinyl  chloride)  above  applicable  MCLs.  Figure  3.4  shows  the  maximum 
concentrations  of  TOC,  dissolved  Mn,  dissolved  Fe,  S2",  and  CH4  at  these  wells;  arsenic  was  not 
monitored  at  any  of  the  wells.  The  concentrations  of  these  SWQIs  were  generally  within 
background  levels.  Where  concentrations  were  much  higher  than  background  (e.g.,  TOC  greater 
than  10  mg/L,  CH4  greater  than  0.1  mg/L),  these  wells  were  generally  located  less  than  10  m 
from  the  injection  zone.  The  production  of  SWQIs  is  therefore  unlikely  to  adversely  impact 
potable  water  supplies. 
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Figure  3.4.  Box  Plots  of  Maximum  TOC,  Dissolved  Mn,  Dissolved  Fe,  Dissolved  S2',  and 
Dissolved  CHj  in  Wells  without  a  Chlorinated  Ethene  Above  Applicable  MCLs. 

Concentrations  of  SWQIs  at  these  wells  were  generally  within  background  ranges,  indicating 
SWQIs  will  not  further  degrade  potable  water  supplies. 
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4.0  SWQI  INDICATOR  SIMULATIONS 

In  this  section,  the  formation  and  natural  attenuation  of  SWQIs  at  ERB  sites  is  illustrated  using  a 
numerical  model  developed  to  simulate  the  transport  and  consumption  of  DOC  and  multiple 
TEAs.  This  model  employs  the  PE  approach  developed  by  Jakobsen  and  Postma  (1999)  and 
demonstrated  by  Curtis  (2003).  This  approach  assumes  that  the  rate-limiting  step  is  fermentation 
of  DOC  (i.e.,  releasing  electron  equivalents)  and  that  other  TEA  reactions  occur  close  to 
equilibrium.  This  model  implementation  is  equivalent  to  organic  carbon  fully  fermenting  to  EE 
and  CO2  at  a  determined  rate,  and  then  H2  becoming  consumed  by  multiple  aqueous  and  mineral 
phase  TEA  processes  controlled  by  their  equilibrium  constants.  As  shown  by  Curtis  (2003),  the 
PE  approach  has  the  advantage  of  avoiding  thermodynamically  unfavorable  reactions  and  large 
errors  in  pH  while  allowing  concomitant  redox  processes  using  multiple  TEAs  (Postma  and 
Jakobsen,  1996;  Jakobsen  and  Postma,  1999;  Jakobsen  et  al.,  1998). 

4.1  MODEL  FORMULATION  AND  INITIAL  CONDITIONS 

The  model  was  used  to  simulate  the  fate  of  the  added  organic  carbon  and  how  the  reduced 
products  affected  the  aquifer  over  a  period  of  40  years  following  addition  of  emulsified  vegetable 
oil  (EVO)  to  form  a  PRB.  The  model  domain  is  a  two-dimensional  cross-section  of  a  water  table 
aquifer  extending  12  m  vertically  and  520  m  horizontally.  The  active  reaction  zone  examined 
with  the  model  extends  260  m  in  the  direction  of  flow.  The  barrier  was  represented  as  a  6  m 
deep  by  15  m  long  by  65  m  wide  zone  containing  a  total  of  4,800  kiloequivalents  (keq)  of 
soybean  oil  (C56H100O6)  evenly  distributed  in  the  treatment  zone  (206  kilograms  per  meter 
barrier  width).  Soybean  oil  dissolution  was  simulated  as  a  first-order  process  so  that  77%  of  the 
residual  soybean  oil  was  released  to  groundwater  as  DOC  in  6.9  years.  This  scenario  results  in 
essentially  all  of  the  soybean  oil  being  released  to  the  aqueous  phase  in  22  years. 

DOC  degrades  by  releasing  electrons  that  react  with  TCE  and  also  with  other  TEAs  including 
O2,  sediment  Mn[IV],  sediment  Fe[III],  SO42",  and  CO2  producing  CH4.  The  reaction  between 
released  electrons  and  TEAs  was  simulated  as  an  equilibrium  reaction.  The  rate  of  DOC 
consumption  (n)  was  simulated  as 

ri  =  -ki  *  Cdoc  -  k2  *  Cdoc  *  Ctce 

where 


Cdoc  =  DOC  concentration  (mM) 

Ctce  =  TCE  concentration  (mM) 
ki  =  1st  order  DOC  decay  rate 
k2  =  2nd  order  DOC  and  TCE  decay  rate 

In  the  partial  equilibrium  framework,  Fe[III]  is  treated  as  an  equilibrium  phase  mineral. 
Dissolution  of  amorphous  Fe(OH)3  and  goethite  (FeOOH)  occur  by  the  reactions: 

Fe(OH)3  +  3H+  Fe3+  +  3H20  pKSp  =  4.9 

and 
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FeOOH  +  3H+  Fe3+  +  2  H20 


pKsp  =  -1.0 


At  many  sites,  Fe[III]  will  be  present  as  mixtures  of  amorphous  Fe(OH)3,  goethite,  and  other 
minerals.  As  a  result,  the  effective  solubility  product  (pKsp)  for  Fe[III]  in  any  aquifer  is 
unknown,  but  is  expected  to  be  somewhere  between  4.9  and  -1.0,  values  representative  of  readily 
dissolved  and  more  stable  Fe  oxide  forms,  respectively.  A  pKsp  value  of  0.239  was  used  in  this 
simulation  based  on  model  calibration  results  from  the  Bemidji  site  and  evidence  from  numerous 
ERB  sites  of  simultaneous  iron  reduction  and  methanogenesis,  which  does  not  occur  when 
higher  dissolution  pKsp  values  are  applied.  Dissolved  Fe[II]  and  Mn[II]  were  assumed  to 
undergo  cation  exchange  with  sediments  and  precipitation  as  siderite  (FeC03)  and  rhodochrosite 
(MnCCL).  Magnetite  precipitation  was  not  considered  in  these  simulations  based  on  limited 
magnetite  production  at  Bemidji. 

CH4  produced  by  fermentation  coupled  to  methanogenic  degradation  initially  enters  the 
dissolved  phase.  Once  the  sum  of  partial  pressures  of  all  dissolved  gasses  (Pch4  +  Pco2  +  Pn2) 
reaches  hydrostatic  pressure,  CH4  and  other  gases  are  assumed  to  outgas  to  the  vadose  zone.  In 
the  field,  outgassing  is  expected  to  occur  in  relatively  coarse  grained  sediments  where  the  gas 
bubbles  can  migrate  upward  through  the  sediment  pores.  In  finer  grained  sediment,  gas  bubbles 
can  remain  trapped  within  the  pores  by  capillary  forces.  Once  CH4  production  slows,  the  CH4  in 
these  trapped  bubbles  could  potentially  dissolve  back  into  groundwater  and  be  carried 
downgradient  with  ambient  groundwater  flow.  In  the  simulations  presented  here,  all  excess  CH4 
is  outgassed  and  there  are  no  trapped  bubbles  that  redissolve. 

The  model  was  run  for  a  total  of  40  years.  The  added  soybean  oil  is  depleted  in  22  years  so  most 
of  the  SWQIs  are  produced  during  the  0-22  year  period  and  then  during  the  22  -  40  year  period 
the  aquifer  plume  continues  to  evolve  and  slowly  recover.  The  initial  water  and  aquifer 
chemistry  for  the  model  is  based  roughly  on  that  of  the  Bemidji  crude  oil  spill  site  (Table  4.1). 
The  model  is  homogeneous  and  isotropic  with  the  groundwater  flow  velocity  set  to  0.06  m/day. 

TABLE  4.1. 


Initial  Condition  and  Background  Concentrations  [mM]  of  Inorganic  Aqueous 
_ Species  and  Sediment  Electron  Acceptors.  _ 


Cone.  (mM) 

Cone.  (mM) 

DIC  (1  mM  =  12  mg/L) 

3.65 

CH4  (1  mM  =  16  mg/L) 

0.0 

Ca2+  (1  mM  Ca  =  40  mg/L) 

1.24 

Cl"  (1  mM  =  35  mg/L) 

0.31 

Fe2+  (1  mM  =  56  mg/L) 

5.1E-21 

Mg2+  (1  mM  =  24  mg/L) 

0.59 

Na+  (1  mM  =  23  mg/L) 

0.077 

02  (1  mM  =  32  mg/L) 

0.25 

Mn2+  (1  mM  =  55  mg/L) 

0 

SO42"  (1  mM  =  96  mg/L) 

1.56 

N2  (inert)  (1  mM  =  28 
mg/L) 

0.45 

PH 

7.6 

Sediment  Fe[III] 

(bioavailable) 

0.98  g/kg 

pe 

14.3 

Cation  Exchange  Capacity 
(CEC) 

47.5  mM/L 
aquifer  volume 

Sediment  Mn[IV] 

0.03  g/kg 
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4.2  SIMULATION  RESULTS 


The  model  results  are  displayed  in  Figures  4.1  and  4.2  as  snapshots  of  the  plume  after  5,  15,  and 
30  years.  The  time  point  at  5  years  was  chosen  to  display  the  SWQI  plume  that  develops  as  a 
result  of  the  ERB  emplacement.  The  15-year  figure  displays  conditions  when  the  soybean  oil  is 
almost  exhausted.  At  this  time,  the  migration  of  the  SWQI  plume  is  apparent,  but  the  waning  of 
the  soybean  oil  source  is  also  clear,  especially  near  the  treatment  zone.  The  final  time  point  of 
30  years  is  used  to  convey  the  progress  of  the  recovery  of  the  aquifer  from  the  reduced 
conditions  created  by  the  remediation. 

Figure  4.1  shows  the  spatial  distribution  of  DOC,  DO,  pH,  SO4,  and  CH4  at  5,  15  and  30  years. 
The  assumed  first  order  dissolution  rate  for  residual  soybean  oil  results  in  65%  being  dissolved 
by  5  years  and  consumption  of  all  the  soybean  oil  after  22  years.  Significant  concentrations  of 
DOC  are  restricted  to  the  area  containing  residual  soybean  oil,  so  CH4  production  is  also 
restricted  to  this  area. 

A  plume  of  depleted  O2  extends  over  100  m  from  the  treatment  zone  by  5  years  and  at  15  years 
extends  beyond  250  m.  Although  recovery  of  the  aquifer  begins  at  22  years,  increased  O2  is 
visible  only  at  the  base  and  top  of  the  model  plume  at  30  years.  The  slow  recovery  of  O 2  is  due 
to  the  large  reservoir  of  sorbed  Fe[II]  on  the  sediments  that  consumes  O2  transported  from 
upgradient.  In  sulfate-rich  aquifers,  FeS  may  form  but  is  not  expected  to  slow  recovery  because 
the  mineral  surface  becomes  oxidized  which  restricts  further  oxidation  of  the  underlying  FeS.  A 
similar  slow  recovery  of  an  aquifer  has  been  observed  at  the  Cape  Cod  research  site  where 
treated  wastewater  entered  the  aquifer  through  infiltration  beds  for  56  years  (Repert  et  al.,  2006). 
While  NO3'  was  not  included  in  the  model  simulation,  NO3"  is  expected  to  behave  similarly  to 
O2  with  slow  recovery  until  the  soybean  oil  and  Fe[II]  are  depleted. 

A  plume  of  slightly  lower  pH  water  migrates  from  the  source  zone.  After  5  years  the  pH  plume 
extends  over  100  m  from  the  source,  and  at  15  years  it  extends  beyond  250  m.  Slightly  higher 
pH  values  can  be  seen  in  the  iron-reducing  zone  due  to  consumption  of  protons  by  the  reaction  of 
DOC  oxidation  coupled  to  iron  reduction.  The  30-year  time  point  shows  that  once  the  soybean 
oil  source  is  exhausted,  pH  in  the  aquifer  recovers  as  clean  groundwater  moves  in  from 
upgradient.  Note  that  the  simulated  pH  changes  are  less  than  one  unit  and  may  not  be  detectable 
because  of  the  typically  sparse  coverage  of  monitoring  wells.  The  pH  is  buffered  by  excess 
calcite  and  cation  exchange  of  H+  with  aquifer  minerals.  More  details  on  the  reactions 
controlling  pH  can  be  found  in  Ng  et  al.  (2015). 

Immediately  after  addition  of  the  soybean  oil,  SO42"  is  depleted  within  the  source  zone.  The 
zone  of  depleted  SO4  extends  approximately  50  m  from  the  source  zone  after  5  years  and 
approximately  75  m  from  the  source  zone  after  15  years.  Fower  SO42"  concentrations  are 
observed  up  to  approximately  140  m  downgradient  after  30  years;  however,  these  levels  are  not 
as  depleted  as  those  closer  to  the  injection  area  in  earlier  time  periods.  In  addition,  after  30  years 
S042"  concentrations  in  and  around  the  source  area  began  to  recover  as  upgradient,  sulfate-rich 
groundwater  fluxes  into  the  area.  Recovery  of  SO42"  is  much  faster  than  O2  since  much  more 
strongly  reducing  conditions  are  required  for  sulfate  reduction. 
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Figure  4.1.  Simulated  Dissolved  Organic  Carbon  (DOC),  Dissolved  Oxygen  (O2),  pH,  Sulfate  (SO4)  and  Dissolved 
Methane  (CH4)  Concentrations  at  5, 15  and  30  Years  after  Substrate  Addition.  DOC,  DO,  SO4,  and  CH4  concentrations 

are  in  units  of  mM. 
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Figure  4.2.  Simulated  Sediment  Fe[III],  Dissolved  Fe2+,  Sediment  Fe[II],  Sediment  Mn[IV],  Dissolved  Mn2+,  and  Sediment 

Mn[II]  at  5, 15  and  30  years  after  Substrate  Addition.  Fe2+  and  Mn2+  concentrations  are  in  units  of  mM.  Sediment  Fe[HI],  Fe[H], 

Mn[IV]  and  Mn[II]  concentrations  are  in  units  of  mmol  per  liter  aquifer  volume. 
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Methane  production  begins  immediately  after  emplacement  of  the  soybean  oil  and  by  5  years  a 
CH4  plume  extends  over  50  m  from  the  source.  By  this  time,  the  total  gas  pressure  in  the 
treatment  zone  has  exceeded  hydrostatic  pressure  and  outgassing  has  started.  By  15  years,  the 
dissolved  CH4  plume  extends  almost  100  m  from  the  source.  The  early  establishment  of 
anaerobic  conditions  in  the  aquifer  after  the  start  of  treatment  restricts  aerobic  CH4  oxidation  to 
small  areas  on  the  fringe  of  the  plume.  Even  after  all  the  soybean  oil  has  been  depleted,  CH4 
consumption  with  oxygen  is  limited  to  a  zone  of  low  O2  and  low  CH4  that  separates  these  two 
plumes. 

In  this  simulation,  CH4  is  assumed  to  be  capable  of  driving  reduction  of  sediment  Fe[III]  based 
on  field  monitoring  results  at  the  Bemidji  site.  If  this  occurs,  then  iron  reduction  will  be  the  main 
process  limiting  the  growth  of  the  dissolved  CH4  plume.  However,  if  adapted  microorganisms 
capable  of  reducing  Fe[III]  using  CH4  are  not  present,  dissolved  CH4  will  migrate  downgradient 
with  little  attenuation  beyond  dispersive  mixing.  Once  the  primary  organic  substrate  (soybean 
oil)  is  exhausted,  CH4  production  stops.  If  essentially  all  of  the  gaseous  CH4  has  been  released 
to  the  vadose  zone,  dissolved  CH4  concentrations  are  expected  to  decline  relatively  quickly. 
However,  if  significant  amounts  of  gaseous  CH4  are  trapped  below  the  water  table,  dissolution  of 
these  trapped  bubbles  could  maintain  the  dissolved  CH4  plume  for  some  time. 

Dissolved  Fe2+  and  Mn2+  and  sediment  Fe[III],  Fe[II],  Mn[IV]  and  Mn[II]  are  shown  in  Figure 
4.2.  Sediment  Fe[II]  and  Mn[II]  includes  both  material  attached  to  cation  exchange  sites  on  the 
sediment  surfaces  and  material  precipitated  as  carbonate  minerals  (siderite  and  rhodochrosite). 

After  5  years  much  of  the  sediment  Fe[III]  is  depleted  in  the  source  zone,  releasing  dissolved 
Fe2+  to  the  aqueous  phase.  However,  downgradient  migration  of  dissolved  Fe2+  is  delayed 
behind  the  low  O2  plume  by  cation  exchange  to  the  sediments.  Dissolved  Fe2+  concentrations 
reach  over  1  mM  while  sediment  Fe[II]  reaches  ~  30  mM  showing  that  most  of  the  reduced  iron 
is  bound  to  the  sediment.  Over  time,  the  simulated  zone  of  depleted  sediment  Fe[III]  grows 
since  the  model  assumes  these  minerals  can  be  reduced  by  dissolved  CH4.  After  all  of  the  DOC 
has  degraded,  the  large  reservoir  of  sorbed  Fe[II]  remains  steady  and  does  not  readily  desorb.  At 
30  years,  the  recovery  of  the  aquifer  is  apparent  as  the  dissolved  Fe[II]  plume  disperses  and 
siderite  redissolves,  but  these  two  forms  make  up  a  relatively  small  fraction  of  the  total  reduced 
Fe.  As  clean  aerobic  groundwater  enters  from  upgradient,  Fe[II]  desorbs  from  the  sediment, 
mixes  with  O2,  and  re-precipitates  as  Fe[III].  During  the  40-year  simulation,  about  25%  of 
bioavailable  Fe[III]  in  the  520-m  model  domain  is  reduced,  but  the  percentage  is  far  higher  in  the 
injection  zone  (-90%). 

Overall,  Mn  behaves  similarly  to  Fe  with  Mn[IV]  reduced  to  dissolved  Mn2+  in  the  source  area, 
and  downgradient  migration  of  Mn2+  limited  by  cation  exchange  on  the  sediments.  The  total 
amount  of  Mn[II]  produced  is  much  lower  than  Fe[II]  due  to  the  much  lower  amount  of  Mn[IV] 
originally  present  in  the  sediment.  This  is  consistent  with  results  from  the  SWQI  database  which 
show  Mn2+  levels  are  often  significantly  lower  than  Fe2+.  The  zone  of  depleted  Mn[IV]  is 
projected  to  grow  more  rapidly  than  the  Fe[III]  depleted  zone  because  Mn  can  be  reduced  under 
relatively  more  oxidizing  conditions  than  Fe.  Similar  to  Fe,  downgradient  migration  of  dissolved 
Mn2+  is  limited  by  cation  exchange  to  the  sediments  and  the  dissolved  Mn2+  dissipates  once  the 
soybean  oil  is  depleted. 
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Figure  4.3  shows  the  distribution  of  electron  (e")  equivalents  released  over  time  in  the  entire 
model  domain  as  soybean  oil  dissolves  releasing  DOC,  which  is  then  fermented  reducing  TCE, 
O2,  Mn,  Fe,  SO42',  and  producing  CH4.  Over  the  first  15  years,  the  very  large  majority  of  e' 
equivalents  are  used  to  reduce  Fe[III]  to  Fe[II],  SO42"  to  S2",  and  CO2  to  CH4.  Overall,  electrons 
consumed  in  O2,  Mn  and  TCE  reduction  are  minor  contributors  to  the  electron  balance.  After  22 
years,  the  fraction  of  e'  equivalents  associated  with  CH4  reduction  declines  slightly  with  a  tiny 
increase  in  Fe[II]  as  CH4  reduces  Fe[III]  in  the  downgradient  aquifer.  As  discussed  above, 
reduction  of  Fe[III]  by  CH4  depends  of  the  presence  of  specialized  microorganisms  that  facilitate 
this  process  and  probably  does  not  occur  at  some  sites. 
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Figure  4.3.  Simulated  Distribution  of  Cumulative  Electron  Equivalents  over  Time 
Following  Soybean  Oil  Addition  to  Stimulate  ERB.  The  contributions  of  the  electron 
acceptors  are  represented  by  the  vertical  heights  of  their  color  bands  at  each  time.  Electron 
equivalents  are  calculated  over  the  first  250  m  downgradient  of  the  injection  point.  The  electron 
equivalents  used  for  TCE  reduction  are  too  low  to  be  visible  on  the  scale  of  the  plot. 


In  summary,  iron  reduction  and  methanogenesis  are  expected  to  be  major  electron  sinks  in  many 
aquifers.  In  aquifers  with  high  levels  of  SO42',  sulfate  reduction  will  also  be  a  major  sink.  Much 
of  the  Fe[II]  and  Mn[II]  produced  is  expected  to  remain  immobilized  on  the  aquifer  material  and 
is  not  expected  to  migrate  long  distances  downgradient.  Similarly,  essentially  all  of  the  sulfide 
will  precipitate  with  available  metals  with  little  to  no  sulfide  migration  downgradient.  If 
appropriate  microorganisms  are  present,  anaerobic  methane  oxidation  using  Fe[HI]  as  an  electron 
acceptor  can  limit  the  downgradient  migration  of  dissolved  CH4.  However  these  organisms  may 
not  be  present,  resulting  in  dissolved  CH4  migrating  substantial  distances  downgradient.  Once 
the  fermentable  organic  carbon  is  depleted,  CH4  production  will  stop  and  the  dissolved  CH4 
plume  is  expected  to  dissipate. 
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5.0  CONCEPTUAL  MODEL  OF  SWQI  PRODUCTION  AND  ATTENUATION 


During  ERB,  large  amounts  of  easily  fermented  organic  substrates  are  added  to  the  target 
treatment  area  to  degrade  or  immobilize  the  COCs.  These  substrates  are  fermented  to  H2, 
acetate,  and  other  VFAs  that  are  then  consumed  in  microbially  mediated  reactions  that  consume 
O2,  NO3",  and  SO42"  as  well  as  Fe[III]  and  Mn[  1 1  I/I  V  ]  containing  minerals,  creating  a  ‘plume’  of 
SWQIs.  Fortunately,  this  ‘plume’  of  impacted  groundwater  is  usually  confined  within  the 
original  contaminant  plume  and  is  unlikely  to  adversely  impact  potable  water  supplies. 

Figure  5.1  shows  a  typical  pattern  of  SWQI  parameters  with  time  in  the:  1)  injection  area;  2) 
near  plume  (25  m  downgradient);  3)  medium-distance  plume  (50  m  downgradient)  and  4)  far 
plume  (100  m  downgradient).  Readers  should  note  that  the  trends  shown  in  Figure  5.1  may  not 
occur  at  all  sites  and  broad  ranges  of  SWQI  parameter  concentrations  have  been  observed.  The 
time  period  for  production  and  attenuation  of  SWQIs  can  vary  from  10  to  over  100  years, 
depending  on  the  amount  and  duration  of  substrate  addition,  groundwater  transport  velocity,  and 
concentrations  of  background  electron  acceptors. 

Organic  substrate  addition  results  in  a  rapid  increase  in  TOC  in  injection  area  monitoring  wells, 
with  maximum  TOC  concentrations  typically  ranging  from  50  to  500  mg/L.  However,  much 
higher  TOC  concentrations  are  observed  at  some  sites.  TOC  concentrations  often  remain  high 
for  several  years  in  the  injection  area  due  to  the  use  of  slow  release  electron  donors  (e.g.,  EVO) 
and/or  repeated  substrate  injections,  and  then  decline  once  substrate  addition  ends.  However, 
low  levels  of  TOC  may  continue  to  be  released  from  endogenous  decay  of  accumulated  biomass 
(Sleep  et  al.,  2005;  Adamson  and  Newell,  2009).  Increases  in  carbon  loading  are  expected  to 
result  in  greater  SWQI  formation.  However,  these  SWQIs  will  attenuate  with  time  and  distance 
downgradient.  Reducing  the  carbon  loading  to  reduce  SWQI  production  may  reduce  treatment 
efficiency,  possibly  resulting  in  greater  exposure  to  chlorinated  solvents  and  other  contaminants. 

Maximum  TOC  concentrations  in  downgradient  wells  are  generally  much  lower  than  in  the 
injection  area,  indicating  TOC  in  the  aqueous  phase  is  rapidly  consumed  and  does  not  migrate 
long  distances  downgradient.  The  rapid  consumption  of  TOC  in  the  injection  area  is  due  to 
reactions  with  background  electron  acceptors  (e.g.  O2,  NO3',  Mn[IV],  Fe[III],  SO42'),  the  target 
contaminants,  and  fermentation  to  CH4.  Since  TOC  is  largely  restricted  to  the  area  in  and 
immediately  downgradient  of  injection  area,  these  redox  reactions  are  also  largely  restricted  to 
the  same  area.  Thermodynamic  calculations  indicate  that  reduction  reactions  should  proceed  in 
the  order  of  O2,  NO3',  Mn[IV],  Fe[III],  SO42',  and  CO2.  However,  these  processes  often  overlap 
(e.g.,  methane  production  occurring  before  complete  sulfate  reduction)  due  to  spatial  variability, 
energy  limitations  from  low  reactant  concentrations,  slow  reaction  kinetics,  and  addition  of 
excess  electron  donor  (Cozzarelli  et  al.,  2000;  Curtis  2003;  Jakobsen  and  Postma  1999). 
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Figure  5.1.  Typical  Variation  in  SWQI  Parameters  over  Time  with  Distance  from 
Injection.  Graphs  compare  concentrations  at  wells  in  the  injection  area  and  25  m,  50  m,  and  100 

m  downgradient  for  40  years  post-injection. 

Declines  in  pH  are  often  observed  during  ERB  due  to  the  release  of  VFAs  and  carbonic  acid 
(H2CO3)  from  substrate  fermentation  and  hydrochloric  acid  (HC1)  from  dechlorination  reactions. 
In  general,  naturally  occurring  carbonate  minerals  are  not  an  effective  pH  buffer,  since  high 
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concentrations  of  dissolved  CO2  produced  during  ERB  prevent  carbonate  mineral  dissolution 
(Kouznetsova  et  al.  2010).  In  some  cases,  pH  declines  can  be  limited  by  alkalinity  released 
during  the  reduction  of  nitrate,  manganese  and  iron  oxyhydroxides,  and  by  CO2  degassing.  If 
sufficient  alkalinity  is  released,  increases  in  pH  may  be  observed. 

In  the  injection  area,  O2  and  NO3'  decline  rapidly  following  substrate  addition  and  often  remain 
low  for  years  after  TOC  declines  due  to  reduction  of  O2  and  NO3  by  sediment  organic  carbon 
and/or  reduced  minerals.  Concentrations  of  O2  and  NO 3'  in  downgradient  wells  decline  with  the 
arrival  of  oxygen-  and  nitrate-depleted  water.  In  most  cases,  there  is  little  or  no  increase  in  O2  or 
NO3'  with  distance  downgradient  due  to  the  very  limited  mixing  between  the  anaerobic  plume 
and  background  aerobic  groundwater  (Borden  and  Bedient  1986;  Schirmer  et  al.  2001). 

Similar  to  O2  and  NO3',  SO42"  concentrations  decline  in  the  injection  area  following  substrate 
addition.  However,  biodegradation  coupled  to  sulfate  reduction  is  less  energetically  favorable 
than  biodegradation  coupled  to  O2,  NO3',  Mn[IV],  and  Fe[III]  reduction.  As  a  result,  SO42'  is 
depleted  more  slowly  than  these  other  TEAs  and  substantial  amounts  of  SO42'  may  persist  in  the 
injection  area  if  background  S042"  levels  are  high.  In  downgradient  wells,  S042"  concentrations 
decline  with  the  arrival  of  anaerobic,  low  SCC2"  groundwater.  However,  the  decline  in  sulfate  in 
downgradient  wells  is  somewhat  more  limited  than  in  the  injection  area,  presumably  due  to 
release  of  SO42'  exchanged  onto  variably  charged  sediment  surfaces  (Bolan  et  al.,  1986; 

Guadalix  and  Pardo,  1991).  Once  injection  area  TOC  declines,  SO42'  is  expected  to  recover 
more  rapidly  than  O2  or  NO3"  since  reduced  iron  and  manganese  minerals  will  not  reduce 
sulfate. 

Sulfate  reduction  results  in  the  production  of  S2-  which  combines  with  protons  (H+)  producing 
hydrogen  sulfide  (H2S).  H2S  is  a  colorless,  poisonous,  inflammable  gas  with  the  characteristic 
‘rotten  egg’  smell  which  is  detectable  at  0.01  ppm  (USEPA  2010).  In  most  aquifers,  substantial 
amounts  of  naturally  occurring  Fe  are  present  and  most  of  the  H2S  produced  will  rapidly 
precipitate  as  insoluble  iron  sulfide  minerals  (FeS  or  FeS2;  Cozzarelli  et  al.  1999).  However, 
detectable  levels  of  sulfide  (~1  mg/L)  are  often  observed  in  injection  areas,  expecially  when 
background  SO42  levels  are  high.  In  the  isolated  cases  where  SO42"  concentrations  are  high  and 
solid  phase  Fe  is  low,  some  dissolved  sulfide  may  migrate  a  short  distance  downgradient  before 
reacting  or  precipitating. 

High  concentrations  of  sulfide  can  inhibit  reductive  dechlorination,  reducing  bioremediation 
effiency.  Hoelen  and  Reinhard  (2004)  reported  that  400  mg/F  Na2S  inhibited  TCE,  DCE,  and 
VC  degradation.  Sung  (2005)  found  that  64  mg/F  total  sulfide  inhibited  growth  of 
de  sulfur  omonas  michiganensis  strain  BB1,  sulfuro spirillum  multivorans,  desulfitobacterium  sp. 
strain  Vietl,  and  dehalococcoides  sp.  strain  FL2  and  strain  BAV1.  However  at  low 
concentrations  (<16  mg/F  as  S),  sulfide  does  not  inhibit  dechlorinators  (Foffler  et  al.,  2005).  At 
most  sites,  sulfide  levels  are  well  below  levels  that  would  inhibit  reductive  dechlorination. 
However,  inhibitory  levels  of  sulfide  were  observed  at  one  site  where  upgradient  sulfate 
concentrations  were  very  high  (-5000  mg/F)  and  bioavailable  iron  was  very  low  (<1 1  mg/kg) 
(Keener  et  al.,  2012). 

Many  aquifer  sediments  contain  significant  levels  of  Fe[III]  and  Mn[  1 1  I/I  V  |  oxyhydroxides. 
Fe[III]  is  commonly  present  as  mixtures  of  poorly  crystalline  amorphous  Fe[III]  oxyhydroxides 
and  more  crystalline  forms.  Solid-phase  Fe  data  compiled  from  20  contaminated  aquifers 
showed  bioavailable  Fe  varied  approximately  3  orders  of  magnitude  from  0.01  g/kg  to  over  17 
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g/kg  and  total  iron  varied  from  0.2  to  1 10  g/kg  (Borden  et  al.,  2015).  Mn  is  most  commonly 
present  as  metastable  Mn[IV]  and  Mn[III]  hydrous  oxides,  characterized  by  structural  defects, 
domain  intergrowths,  cation  vacancies,  and  solid  solutions  (Villalobos  et  al.  2003).  The 
Mn[III/IV]  concentration  in  soils  and  aquifer  sediments  (median  =  0.02  g/kg;  Chen  et  al.,  1999) 
is  often  much  lower  than  the  iron  concentration. 

During  ERB,  increased  levels  of  Eb  and  VFAs  produced  from  substrate  fermentation  lead  to 
reduction  of  Fe[III]  and  Mn[  1 1  I/I  V  ]  oxyhydroxides  and  release  of  dissolved  Fe2+  and  Mn2+. 

Early  research  suggested  that  the  poorly  crystalline  forms  of  Fe[III]  (e.g.,  ferrihydrite)  could  be 
rapidly  and  completed  reduced  (Lovley  and  Phillips  1986)  while  reduction  of  more  crystalline 
forms  (e.g.,  goethite  and  hematite)  was  more  limited  (Lovley  and  Phillips,  1987).  However, 
more  recent  work  has  shown  that  a  variety  of  crystalline  Fe-oxides  can  be  reduced  by  ERB 
(Roden  and  Zachara,  1996;  Roden  et  al.,  2000;  Roden  and  Urrutia,  2002).  Mn[II]  can  be 
generated  by  dissimilatory  reduction  of  Mn(IIEIV)  hydrous  oxides  on  sediments  by 
microorganisms  (Tebo  et  al.,  2004)  or  by  abiotic  reduction  of  these  oxides  by  Fe[II]  (Postma  and 
Appelo,  2005).  In  some  aquifers,  Mn  and  Fe  concentrations  have  been  observed  to  decline 
before  TOC  is  depleted,  presumably  due  to  depletion  of  bioavailable  Mn[  III/IV]  and  Fe[III]  in 
the  injection  area  aquifer  material.  Once  TOC  levels  decline  in  the  injection  area,  production  of 
additional  Mn[II]  and  Fe[II]  will  slow,  and  dissolved  Mn  and  Fe  concentrations  in  monitoring 
wells  should  start  to  decline. 

Much  of  the  reduced  Fe[II]  produced  during  ERB  will  be  retained  on  the  aquifer  material  within 
the  injection  area  through  a  variety  of  processes  including  ion  exchange  and  surface 
complexation  reactions  (Davis  and  Kent,  1990;  Baeyens  and  Bradbury,  1997;  Davis  et  al.  1998) 
and  formation  of  mixed  Fc[  II/IIIJ  phases  such  as  magnetite  and  green  rust  and  pure  Fe[II]  phases 
like  siderite  (Frederickson  et  al.,  1998;  Hansel  et  al.  2005).  Mn[II]  is  not  expected  to  precipitate 
as  most  sites,  but  can  be  retained  by  sorption  processes  when  concentrations  of  competing 
cations  are  low  and  the  pH  is  alkaline  (Baeyens  and  Bradbury,  1997).  However,  ERB  often 
results  in  reduced  pH  with  high  concentrations  of  Fe2+  which  will  minimize  Mn2+  retention  in  the 
injection  area  aquifer  material. 

Fe[II]  and  Mn[II]  that  is  not  retained  within  the  injection  area  will  migrate  downgradient  with 
natural  groundwater  flow.  During  downgradient  migration,  a  portion  of  the  dissolved  Mn  and  Fe 
will  sorb  to  the  aquifer  material,  reducing  the  maximum  concentrations  observed  in 
downgradient  wells.  Over  time,  Fe2+  and  Mn2+  production  in  the  source  area  will  decline,  and 
these  sorbed  materials  could  potentially  be  released  to  the  aqueous  phase.  However,  the 
desorption  rate  will  be  limited  by  the  supply  of  exchangeable  cations  entering  the  plume  from 
upgradient.  Numerical  model  simulations  of  oxygen  recovery  in  an  anaerobic  wastewater  plume 
in  Cape  Cod,  MA  indicate  that  the  maximum  concentrations  of  dissolved  Fe  and  Mn  released 
will  be  much  lower  than  peak  concentrations  observed  during  the  substrate  loading  phase  (Ng 
2015). 

Arsenic  is  naturally  present  in  many  aquifers,  often  co-precipitated  and/or  sorbed  to  iron  and 
aluminosilicate  minerals  (Pierce  and  Moore,  1982;  Smedley  and  Kinniburgh,  2002;  Kent  and 
Fox,  2004).  The  general  term  sorption  is  used  to  encompass  all  binding  mechanisms  of  arsenic 
to  Fe[III]  phases.  Under  oxidizing  conditions,  As[V]  is  the  dominant  form  and  is  present  in 
solution  as  a  negatively  charged  ion  (HoAsOT  or  HAsCF2").  Under  mildly  reducing  conditions, 
As[III]  is  the  dominant  form  and  occurs  as  an  uncharged  species  (As(OH)3°).  Typically,  As[III] 
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sorbs  less  strongly  than  As[V]  on  minerals  except  Fe[III]  oxides  and  oxyhydroxides  (Dzombak 
and  Morel,  1990;  Dixit  and  Hering,  2003). 

During  ERB,  sorbed  arsenic  is  released  to  solution  when  As[V]  is  reduced  to  As  [III]  and  Fe[III] 
is  reductively  dissolved  (Kent  and  Fox,  2004;  Hohn  et  al.,  2006;  Hering  et  al.,  2009;  He  et  al., 
2010;  Cozzarelli  et  al.,  2015).  As  Fe-rich  groundwater  migrates  downgradient,  Fe[II]  sorbs  to 
the  sediment  or  precipitates  as  Fe-bearing  minerals  (FeS,  carbonates,  magnetite)  and  aqueous  Fe 
concentrations  decline.  The  available  monitoring  data  suggest  that  arsenic  follows  a  similar 
pattern  and  aqueous  arsenic  concentrations  decline  as  the  anaerobic  plume  migrates 
downgradient  and  encounters  Fe[III]-rich  sediments.  While  oxidation  of  aqueous  As[III]  to 
As[V]  is  a  kinetically  slow  process,  As[III]  can  sorb  to  amorphous  ferrihydrite  (am-Fe(OH)3) 
and  then  undergo  much  more  rapid  heterogeneous  oxidation  (i.e.,  oxidation  at  an  iron-mineral 
surface),  generating  strongly  sorbed  As[V]  (Pierce  and  Moore,  1982;  Smedley  and  Kinniburgh, 
2002).  Once  TOC  concentrations  in  the  injection  area  decline,  Fe[III]  and  As[V]  reduction  is 
expected  to  decline  with  a  concurrent  decline  in  arsenic  release.  At  sites  where  concentrations  of 
sediment-bound  As[V]  are  low,  minimal  arsenic  will  be  released. 

TOC  fermentation  products  (H2  and  acetate)  in  the  injection  area  that  are  not  consumed  in 
microbially  mediated  reactions  with  contaminants  or  background  electron  acceptors  will  be 
fermented  to  CH4,  which  can  result  in  high  CH4  concentrations  in  the  injection  area  and 
downgradient  aquifer  (Jacob  et  al.,  2005;  Riis  et  al.,  2007).  In  saturated  media,  if  the  sum  of  gas 
partial  pressures  exceeds  the  hydrodynamic  pressure,  the  gases  will  come  out  of  solution  and 
form  bubbles  (Amos  and  Mayer,  2006).  This  occurs  primarily  near  the  water  table  due  to  the 
lower  hydrostatic  pressure.  At  greater  depths,  groundwater  can  appear  to  be  supersaturated  with 
CH4  due  to  the  higher  pressure.  In  relatively  homogeneous,  coarse  grained  sediments,  gas 
bubbles  can  migrate  upward  into  the  vadose  zone,  removing  CH4  from  the  aquifer.  However,  in 
finer  grained  sediments  upward  migration  of  gas  bubbles  will  be  more  limited.  If  the  gas 
bubbles  are  not  released  to  the  vadose  zone,  they  will  eventually  dissolve,  releasing  CH4  back 
into  groundwater. 

Field  monitoring  at  ERB  sites  has  shown  that  CH4  levels  can  be  elevated  in  the  headspace  of 
injection  area  wells,  while  CH4  levels  in  soil  gas  measurements  collected  just  below  the  land 
surface  are  typically  lower  (TSI  and  Solutions-IES,  2003).  This  is  due  to  diffusive  mixing  with 
O2  from  the  land  surface  and  subsequent  oxidation  by  methanotrophs  (Fischer  et  al.,  1996;  Amos 
et  al.,  2005;  Ma  et  al.,  2012;  Warren  et  al.,  2014).  However,  CH4  can  be  released  to  the  land 
surface  or  subsurface  structures  when  methane  production  rates  are  high  and  oxygen  transport 
through  the  vadose  zone  is  restricted  by  saturated  soil  conditions  (Sihota  et  al.,  2013). 

Dissolved  CH4  produced  in  the  injection  area  will  migrate  downgradient  with  groundwater  flow 
and  could  potentially  biodegrade  if  the  CH4  plume  mixes  with  aerobic  groundwater.  However, 
the  actual  extent  of  aerobic  CH4  biodegradation  in  the  saturated  zone  will  be  limited  by  an 
anoxic  zone  that  develops  between  the  strongly  anaerobic  water  containing  CH4  and  the 
oxygenated  background  groundwater  (Borden  and  Bedient  1986,  Amos  et  al.  2005).  Recent 
research  indicates  that  methane  could  be  consumed  using  NO 3'  (Raghoebarsing  et  al.,  2006),  Mn 
and  Fe  oxides  (Beal  et  al.,  2009;  Crowe  et  al.,  2011;  Amos  et  al.  2012)  and  SO42"  (Caldwell  et 
al.,  2008;  Grossman  et  al.,  2002)  as  TEAs.  These  anaerobic  TEAs  are  potentially  very  important 
since  they  could  limit  the  downgradient  migration  of  CH4  plumes.  However  anaerobic  methane 
oxidation  requires  a  consortium  of  archaea  and  bacteria  (Boetius  et  al.,  2000)  and  has  a  very  low 
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energy  yield  (Caldwell  et  al.,  2008).  As  a  result,  anaerobic  methane  oxidation  reactions  may  be 
slow  or  may  not  occur  at  some  sites,  and  dissolved  CH4  can  migrate  long  distances  in  some 
aquifers.  Once  TOC  in  the  injection  area  declines,  CH4  production  will  stop  and  dissolved  CH4 
should  be  transported  downgradient  by  groundwater  flow. 
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6.0  RESEARCH  IMPLICATIONS  AND  OPPORTUNITIES 


At  most  sites,  the  majority  of  the  ED  added  to  the  aquifer  is  consumed  in  the  reduction  of  Fe[III] 
to  Fe[II],  S042‘  to  S2",  and  CO2  to  CH4.  At  many  sites,  anaerobic  methane  oxidation  using 
Fe[III]  and  SO4  as  TEAs  is  expected  to  be  the  primary  attenuation  mechanism  for  dissolved 
CH4,  resulting  in  further  production  of  Fe[II]  and  S2‘  at  many  sites.  The  very  large  majority  of 
Fe[II]  and  S2'  produced  is  retained  on  the  aquifer  material,  generating  a  large  reservoir  of 
reducing  power  within  the  aquifer. 

Monitoring  data  and  modeling  studies  at  the  Cape  Cod,  MA  wastewater  plume  indicate  that  a 
portion  of  the  added  organic  carbon  is  retained  on  the  aquifer  material  as  BMSOC.  The  BMSOC 
is  a  small  fraction  of  sediment  organic  carbon  and  oxidation  of  BMSOC  appears  to  be  slow. 
However,  this  BMSOC  may  greatly  delay  re-oxidation  of  the  aquifer  and  has  the  potential  to 
generate  reducing  conditions  that  support  further  Fe[II]  reduction. 

These  large  reservoirs  of  reduced  iron  and  sulfur  minerals  and  organic  carbon  will  cause  the 
aquifer  downgradient  of  ERB  projects  to  remain  anoxic  to  reducing  for  years  to  decades  after  the 
end  of  substrate  addition.  However,  these  reservoirs  of  Fe[II],  reduced  S,  and  organic  carbon 
also  have  the  potential  to  increase  degradation  and/or  attenuation  rates  for  a  variety  of 
groundwater  contaminants  including  chlorinated  solvents,  explosives  and  propellants  (TNT, 
RDX,  HMX,  C104").  NO3',  chromate,  and  selected  radionuclides  (Tc04',  U022+)  (Butler  and 
Hayes,  1998;  Lee  and  Batchelor,  2002,  2003;  Ferry  et  al.,  2004;  Sleep  et  al.,  2005;  AFCEE, 

2008;  He  et  al.,  2009;  Adamson  and  Newell,  2009).  Additional  research  is  needed  to  better 
understand  the  impact  of  these  reduced  materials  on  contaminant  degradation/attenuation  rates. 
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Table  A-l.  Characteristics  of  ERB  Sites  Included  in  SWQI  Database, 


No.  of  Record 


Site  ID 

Location 

Injection  Type 

Organic  Substrate 

Data  Source 

wells 

(months) 

Aerojet  General  Corporation 

Rancho  Cordova,  CA 

Horizontal  Wells 

Citric  Acid 

ESTCP  Project  Report 

18 

24 

Aiea  Laundry 

Pearl  Harbor,  HI 

Source 

EOS,  Lactate 

Personal  Correspondence 

2 

6 

Altus  AFB  OU1  -  Biowall 

Altus,  OK 

Biowall 

Mulch  and  compost 

AFCEE  Project  Report 

44 

88 

Altus  AFB  OU1  -  Injections 

Altus,  OK 

Source 

EVO 

Monitoring  Report 

5 

20 

Altus  AFB  SS-17 

Altus,  OK 

PRB 

EVO 

Monitoring  Report 

14 

20 

Altus  AFB  SS-18 

Altus,  OK 

Source 

EOS 

Monitoring  Report 

22 

21 

Altus  AFB  SS-22 

Altus,  OK 

Source 

EOS 

Monitoring  Report 

52 

21 

Avon  Park  AFB 

Avon  Park,  FL 

Recirc 

Lactate,  ethanol 

SERDP  Project  Report 

5 

5 

Beale  AFB  Site  39  Building  2145  Marysville,  CA 

Source 

EOS 

Monitoring  Report 

9 

15 

Beale  AFB  Site  39  Source  Area 

1 

Marysville,  CA 

Source 

EOS 

Monitoring  Report 

18 

10 

Boeing  Former  Compton  Site 

Compton,  CA 

Source 

Newman  Zone 

Monitoring  Report 

14 

36 

Cape  Canaveral  Flangar  K 

Cape  Canaveral,  FL 

Source 

Soybean  oil 

AFCEE  Project  Report 

19 

68 

Charleston  NWS  SWMU  17 

Charleston,  SC 

Source 

EOS 

ESTCP  Project  Report 

26 

41 

Cornhusker  AAP  OU1 

Grand  Island,  NE 

PRB 

Newman  Zone,  Wesblend  66 

Monitoring  Report 

47 

48 

Crown  Cork  &  Seal  Facility 

Pittsburg,  CA 

Source 

HRC 

Monitoring  Report 

2 

66 

Dover  AFB  Area  5 

Dover,  DE 

PRB 

Lactate,  Newman  Zone 

Monitoring  Report 

43 

77 

Dover  AFB  Area  6 

Dover,  DE 

PRB 

Lactate,  Newman  Zone 

Monitoring  Report 

136 

77 

Dover  AFB  Site  WP14 

Dover,  DE 

Biowall 

Mulch 

AFCEE  Project  Report 

14 

39 

East  Charleston  Business  Park 

Mountain  View,  CA 

PRB 

HRC 

Monitoring  Report 

9 

64 

Fort  Devens  AOC  50 

Devens,  MA 

PRB 

Molasses 

Monitoring  Report 

9 

30 

Fort  Dix  MAG-1 

Fort  Dix,  NJ 

Recirc 

Lactate 

ESTCP  Project  Report 

22 

12 

Fort  Worth  Joint  Reserve  Base 
AOC-2 

Fort  Worth,  TX 

PRB 

Soybean  oil  and  HF  corn  syrup 

AFCEE  Project  Report 

14 

26 

Plickam  AFB  Site  LF05  - 

Southeast  Area 

Honolulu,  HI 

Recirc 

Lactic  acid 

Monitoring  Report 

37 

4 

Plickam  AFB  Site  LF05  -  TCRA 

Area 

Honolulu,  HI 

PRB 

Soybean  oil 

AFCEE  Project  Report 

11 

26 

Indian  Head  Naval  Surface 
Warfare  Center 

Indian  Head,  MD 

Recirc 

Lactic  acid 

Literature  Report 

14 

5 

JEB  Little  Creek  Site  1 1 

Virginia  Beach,  V A 

PRB 

SRS  emulsified  oil 

Monitoring  Report 

26 

12 

JEB  Little  Creek  Site  12 

Virginia  Beach,  VA 

PRB 

EOS 

Monitoring  Report 

17 

13 

JEB  Little  Creek  Site  13 

Virginia  Beach,  VA 

PRB 

LactOil,  ethyl  lactate 

Monitoring  Report 

28 

6 
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Table  A-l  (continued).  Characteristics  of  ERB  Sites  Included  in  SWQI  Database. 


Data 

No.  of  Record 


Site  ID 

Location 

Injection  Type 

Organic  Substrate 

Data  Source 

wells 

(months) 

Maryland  Perchlorate  Site 

Elkton,  MD 

PRB 

EOS 

ESTCP  Project  Report 

8 

42 

McConnell  AFB  Site  FT07 

Wichita,  KS 

Biowall/Recirc 

Mulch  and  soybean  oil 

Monitoring  Report 

18 

8 

Mercedes  North  Houston 

Houston,  TX 

Source 

HRC,  Newman  Zone,  Cheese  Whey 

Monitoring  Report 

10 

51 

Moffett  Field  Site  26 

Mountain  View,  CA 

Source 

EHC 

Monitoring  Report 

14 

22 

Moffett  Field  Site  28  -  Building 

88 

Mountain  View,  CA 

Source 

Lactate  and  lactic  acid 

Monitoring  Report 

7 

10 

Moffett  Field  Site  28  -  Traffic 
Island 

Mountain  View,  CA 

Source 

Lactoil 

Monitoring  Report 

6 

10 

Moffett  Field  Site  28  -  W9-18 
Area 

Mountain  View,  CA 

Source 

EHC 

Monitoring  Report 

12 

10 

Naval  Weapons  Station  Seal 
Beach  -  Site  40 

Seal  Beach,  CA 

Source 

Sodium  lactate,  HRC 

Monitoring  Report 

7 

75 

Oates  Park 

Mesquite,  TX 

Source 

HRC 

Monitoring  Report 

34 

36 

Pacific  Scientific 

Santa  Barbara,  CA 

Source 

HRC 

Monitoring  Report 

4 

110 

Picatinny  Arsenal  Area  157 

Dover,  NJ 

Recirc 

Cheese  whey 

ESTCP  Project  Report 

14 

18 

Renco  Encoders 

Goleta,  CA 

Source,  PRB 

HRC,  HRC-X,  Wilclear,  Lactoil,  EOS,  EHC 

Monitoring  Report 

14 

117 

Savannah  River  Site  -  F  Area 

Aiken,  SC 

Source 

Molasses 

Monitoring  Report 

27 

19 

Tinker  AFB  FTA-2 

Oklahoma  City,  OK 

Source 

Soybean  oil  and  HF  corn  syrup 

AFCEE  Project  Report 

20 

23 

Travis  AFB  Site  DP039 

Fairfield,  CA 

Biowall/Recirc 

Mulch  and  soybean  oil 

Monitoring  Report 

10 

11 

Treasure  Island  Site  24 

San  Francisco,  CA 

Recirc 

Sodium  lactate 

Monitoring  Report 

9 

12 

TRW  Microwave 

Santa  Clara,  CA 

Source,  PRB 

HRC,  Cheese  Whey,  EVO,  EHC-L,  ABC+ 

Monitoring  Report 

20 

132 

Western  Microwave 

Sunnyvale,  CA 

PRB 

HRC 

Monitoring  Report 

33 

68 

Whittaker  Bermite  Facility  Area 

Santa  Clarita,  CA 

Recirc 

Citric  Acid 

Monitoring  Report 

13 

7 

Notes 

PRB  -  Permeable  reactive  barrier 
Recirc  -  Recirculation  system 
Source  -  Source  area  injections 
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Table  A-2.  Summary  Statistics  for  Post-Injection  Data  Contained  in  SWQI  Database 


Well  Location 

Monitoring 

Percentiles 

Points 

Min 

10th 

25th 

50th  | 

75th  | 

90th 

Max 

Maximum 

Total 

Organic 

Carbon 

(mg/L) 

Upgradient 

113 

<1 

1.4 

2.3 

3.8 

12  j 

38 

258 

Injection  Area 

303 

0.3 

3.1 

12 

230  | 

1,038  i 

3,548 

33,000 

0-10  m  Downgradient 

174 

<1 

3.2 

10 

57  i 

263  j 

738 

7,890 

10-25  m  Downgradient 

64 

<1 

1.9 

2.8 

7.0  ! 

81  j 

385 

5,080 

25-50  m  Downgradient 

35 

<1 

1.6 

<5 

9.6  ! 

20  i 

125 

885 

50+  m  Downgradient 

75 

0.833 

2  2 

3.0 

5.7  | 

15  1 

28 

864 

Upgradient 

61 

0.002 

0.0066 

0  057 

0.2  i 

0.53  I 

1.1 

39 

Maximum 

Injection  Area 

125 

0.023 

0.15 

0.80 

2.6 

8.7  j 

17 

53 

Dissolved 

0-10  m  Downgradient 

112 

0.0065 

0.018 

0.087 

1.1  i 

4.5  ; 

12 

86 

Manganese 

10-25  m  Downgradient 

52 

0.0015 

0.020 

0.13 

0  44  j 

1.7 

7.1 

59 

(mg/L) 

25-50  m  Downgradient 

20 

0.002 

0  010 

0.024 

0.077  1 

0.20  ! 

0.97 

2.0 

50+  m  Downgradient 

17 

0.006 

0  018 

0042 

0.2  i 

3.2  | 

15 

31 

Upgradient 

109 

<0.04 

<0.1 

0.1 

0.43 

1.1  i 

6.0 

350 

Maximum 

Dissolved 
Iron  (mg/L) 

Injection  Area 

203 

0.02 

0.14 

2.0 

22  1 

120  | 

266 

1,200 

0-10  m  Downgradient 

180 

0.02 

0.16 

1.2 

8.2  j 

78  i 

251 

928 

10-25  m  Downgradient 

64 

<0.1 

0.17 

0.50 

6.9  ! 

50  j 

218 

470 

25-50  m  Downgradient 

33 

<0.1 

0.38 

1 

8.9  j 

75 

167 

404 

50+  m  Downgradient 

69 

0.01 

0.1 

0.2 

0.74  ! 

5  1 

21 

200 

Upgradient 

39 

<0  00031 

0.00031 

0.0016 

0.0035  i 

0.005  ! 

0.012 

0.073 

Maximum 

Injection  Area 

44 

0.0011 

<0.01 

0.0088 

0.033  j 

0  11  j 

0.49 

1.2 

Dissolved 

0-10  m  Downgradient 

60 

0.0011 

0.0028 

<0.01 

0.014  ! 

0.040  | 

0.096 

0.90 

Arsenic 

10-25  m  Downgradient 

18 

0  0012 

0.0040 

<0.012 

0.011  | 

0.042  I 

0.36 

0.75 

(mg/L) 

25-50  m  Downgradient 

6 

<0.005 

<0  0075 

<0.01 

0.0056  i 

0  022 

0.029 

0.032 

50+  m  Downgradient 

10 

<0.005 

<0.005 

<0.005 

0.0056  i 

0.06768  j 

0.111 

0.12 

Upgradient 

124 

<0.1 

4.5 

14 

81 

281  i 

1,285 

1,870 

Minimum 

Sulfate 

(mg/L) 

Injection  Area 

300 

<0.1 

0.36 

1 

9.7  | 

81 

450 

3,150 

0-10  m  Downgradient 

222 

<0.2 

<0.5 

<0.2 

8.5  I 

98  i 

843 

3,100 

10-25  m  Downgradient 

78 

0.12 

0.49 

1.5 

13 

96  j 

643 

2,930 

25-50  m  Downgradient 

36 

0.1 

<0.5 

0.70 

11  1 

120  | 

210 

2,900 

50+  m  Downgradient 

75 

0.1 

0.37 

1 

13  j 

120  i 

416 

2,400 

Upgradient 

59 

0.0215 

<0.05 

0.0745 

<1  j 

0.575  ! 

2.32 

190 

Maximum 

Sulfide 

(mg/L) 

Injection  Area 

119 

<0.05 

0.256 

0.78 

2  j 

10 

110.2 

344 

0-10  m  Downgradient 

86 

0.01 

0.086 

0.39375 

1.06  i 

2.21  i 

9.55 

1300 

10-25  m  Downgradient 

38 

<0.04 

<0.1 

0  4025 

0.88  j 

2  535  : 

8.144 

20 

25-50  m  Downgradient 

21 

<0.05 

<0.1 

0.3 

<1 

0.9  j 

1.3 

1.6 

50+  m  Downgradient 

15 

<0.1 

<1 

0.6 

<2  j 

1.75  ! 

2.54 

2.9 

Upgradient 

107 

0  0001 

0  0015 

0.0029 

0.063 

0.32  j 

1.3 

6.2 

Maximum 

Methane 

(mg/L) 

Injection  Area 

261 

<0  00031 

0.050 

1.2 

9.0  j 

15  i 

23 

120 

0-10  m  Downgradient 

197 

0.0028 

0.20 

0.91 

7.6  j 

17 

31 

77 

10-25  m  Downgradient 

69 

0.00037 

0.0056 

<0.025 

1.9 

11  ! 

28 

56 

25-50  m  Downgradient 

36 

0.0002 

0.0045 

0.059 

1.9  ! 

17  i 

24 

40 

50+  m  Downgradient 

74 

0.00022 

<0.002 

0.013 

0.25 

4.9 

17 

35 

51 


Table  A-2  (continued).  Summary  Statistics  for  Post-Injection  Data  Contained  in  SWQI  Database, 


Well  Location 

Monitoring 

Percentiles 

Points 

Min 

10th 

25th 

;  50th  ; 

75th 

90th 

Max 

Upgradient 

130 

0 

0 

0.02 

1  0.32  I 

1.4 

2.9 

7.3 

Minimum 

Injection  Area 

323 

0 

0 

0.05 

j  0.15  i 

0.45 

1.1 

8.8 

Dissolved 

0-10  m  Downgradient 

208 

0 

0 

0 

1  0.07  ! 

0.2 

1.3 

10 

Oxygen 

10-25  m  Downgradient 

81 

0 

0 

0.04 

I  0  18  ! 

0.81 

1.8 

6.1 

(mg/L) 

25-50  m  Downgradient 

39 

0 

0 

0.01 

j  0.11  i 

0.28 

0.60 

1.4 

50+  m  Downgradient 

79 

0 

0.07 

0.12 

!  0.2  | 

0.46 

1.6 

9.4 

Upgradient 

73 

0.0011 

0.028 

<0.1 

<0.22  ! 

1.7 

4.8 

12 

Minimum 
Nitrate  (mg/L 
as  N) 

Injection  Area 

136 

<0.0066 

00058 

<0  05 

I  <0.11  j 

0.19 

0.59 

14 

0-10  m  Downgradient 

132 

0  00068 

00055 

<0  05 

1  <0.2  j 

0.20 

0  41 

3.5 

10-25  m  Downgradient 

45 

<0.01 

<0.01 

<0  05 

|  <02  | 

0.90 

6.2 

25 

25-50  m  Downgradient 

16 

<0.01 

<0.01 

<001 

|  0.075  j 

0.51 

2.2 

2.5 

50+  m  Downgradient 

21 

<0.01 

<0.05 

<0  11 

j  o  16  j 

1.3 

2.5 

5.9 

Upgradient 

106 

3.9 

5.8 

6.6 

1  6.8  ; 

7.1 

7.2 

11.7 

Median  pH 

Injection  Area 

223 

4.9 

5.7 

6.2 

6.6  ! 

6.8 

72 

8.7 

(Non- 

0-10  m  Downgradient 

132 

4.8 

6.1 

6.4 

6.7 

7.0 

7.2 

8.4 

Buffered) 

10-25  m  Downgradient 

55 

4.9 

6.3 

6.7 

j  7.0  j 

7.2 

7.5 

10.0 

(SU) 

25-50  m  Downgradient 

16 

6.1 

6.6 

6.8 

j  7.0  j 

7.2 

7.6 

7.7 

50+  m  Downgradient 

46 

5.8 

6.2 

6.6 

6.9 

7.1 

7.2 

9.4 

Upgradient 

41 

3.6 

4.6 

4.8 

j  5.4  1 

5.9 

6.6 

7.2 

Median  pH 
(Buffered) 
(SU) 

Injection  Area 

106 

4.2 

4.9 

5.3 

!  6.0  | 

6.5 

6.8 

8.9 

0-10  m  Downgradient 

86 

3.5 

5.4 

5.8 

i  6.1  i 

6.5 

6.9 

10.3 

10-25  m  Downgradient 

25 

3.8 

5.1 

5.9 

j  6.1  i 

6.5 

6.8 

11.2 

25-50  m  Downgradient 

23 

4.7 

5.1 

5.3 

5.8  ; 

6.2 

6.5 

7.7 

50+  m  Downgradient 

32 

3.7 

4.8 

4.9 

!  5.1  ; 

5.4 

6.4 

8.2 

Upgradient 

76 

-1.9 

-0.90 

-0.32 

i  -0.12  i 

0.02 

0.25 

1.2 

Median  ApH 

Injection  Area 

165 

-3.3 

-1.6 

-0.90 

i  -0.35  i 

000 

0.31 

1.4 

(Non- 

0-10  m  Downgradient 

108 

-3.8 

-1.2 

-0.54 

!  -021  i 

0.15 

0.63 

1.0 

Buffered) 

10-25  m  Downgradient 

47 

-3.2 

-1.2 

-0.58 

!  -0.20  ! 

0.01 

0.19 

2.7 

(SU.) 

25-50  m  Downgradient 

9 

-0.68 

-0.49 

-0.39 

1  -0.21  I 

0.00 

0.11 

0.38 

50+  m  Downgradient 

34 

-1.6 

-0.33 

-0.17 

-0.10  j 

0.02 

0.07 

0.24 

Upgradient 

35 

-0.90 

-0.42 

-0.20 

I  o.io  ; 

0.44 

081 

1.4 

Median  ApH 
(Buffered) 
(SU.) 

Injection  Area 

57 

-2.0 

-1.3 

-0.4 

!  04  j 

1.3 

2.1 

4.4 

0-10  m  Downgradient 

73 

-2.9 

-0.13 

0.19 

|  0.47  1 

1.0 

1.6 

3.7 

10-25  m  Downgradient 

16 

-2.8 

-1.2 
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